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ABSTRACT 


This thesis presents some general equations describing the 
particle density distributions within a plasma. Solutions to 
these equations are discussed for various geometries, with particular 
emphasis given to the effects of diffusion losses at probes of dimen- 
sions comparable to or greater than the mean free path and to magnetic 
field end effects. For many of the geometries discussed, including 
one probe geometry, the effects of a superimposed magnetic field is 


discussed, 





INTRODUCT ION 


One of the complex problems in the field of plasma physics is 
that of predicting the effects of surfaces with which the plasma is in 
contact upon the particle density distribution functions and related 
phenomena. No plasma is infinite, and consequently a plasma must be 
confined by walls, which may be either conductors or insulators. 
Further, in order to make measurements within a plasma, probes may 
be inserted. What is their effect on the particle density distri- 
bution and related phenomena such as potential distribution and 
energy dissipation? What are the particle currents to various surfaces? 
The latter are particularly important because they determine the 
heat liberated at the surface, which in turn is an energy loss which 
may cause very high temperatures, increased local ionization, and 
melting of materials in extreme cases. This thesis discusses these 
effects when the surfaces are insulators. 

In this thesis some general equations describing a plasma are 
presented with their corresponding restrictions. Solutions to 
these equations are discussed for various geometries, with particular 
emphasis given to the effects of probes and to magnetic field end 
effects on the steady-state density distributions and on the energy 
losses at the surfaces involved. For many of the geometries dis- 
cussed, including one probe geometry, the effects of a superimposed 


magnetic field is studied. 





ee 


GENERAL RELATIONS 
In this section, general equations describing a plasma are 
developed and discussed. 


Ae Conservation Equations 


The Boltzmann Transport Equation 


Zl) HOR) = mE WK - nx TIA aff FO pARVArr 


is used as the starting point for the development of the macroscopic 
relations governing the plasma. By ascribing appropriate values to X 


and x* {1} (2], we obtain the continuity equation: 


ao +Ve(nyy) = nO -y) Ce) 
(in which y, is the mean ionizing frequency and )), the mean recombination 


frequency) , and the equation for conservation or momentum: 


~ A V2 as Sa - —_ oh 
Se(nv) + i —n(&+OxVa) = -nz V4 (2=3) 


in which the right hand term is really a definition of the mean 


—- 


collision frequency, VY, , the pressure ‘Dp is a symmetric tensor 


Pi = NM V- » ( 2=3a) 
and in which 
Dae ae (2=3b) 
mM 
and 
lu, = — * e (2eee) 


_> = 
E, By qy my and vq are, respectively, the electric field intensity, 
the magnetic field intensity, particle charge, particle mass, and 


Griit Velocity. 





The continuity equation (2-2) can be further simplified. We 
shall consider a two component system of electrons and singly=charged 
ions in which volume ionization by electrons is essential. Recombina- 
tion takes place predominately at the walls, with volume recombination 


being negligible. We have, then, 


Ons OT 
ar + V- J". =—eait= (2-);) 





where Y is the net ionizing frequency and the particle current density 


is 


(aes, ( 2-la) 
In the above equations, and in the remainder of this thesis, the 
subscript "+" refers to singly=charged positive ions, and the sub- 
script "=!" refers to electrons. Under those conditions for which 
the rates of decay of ions and electrons are eSsentially equal, so 
that 
i = (2-5) 
we obtain for the continuity equaticn 
Val. <a =n. - 2 (2-6) 
In the momentum equation (2=3), the pressure term is generally 
a difficult one to evaluate because ‘p is not known, and the 
cependence on its divergence makes good approximations difficult. 
However, for the simple case in which the distribution of random 


velocities ts essentially isotropic, so that 


VV & O On. (2=7a) 
and 
V&V, 8s, (2=75) 





the pressure tensor reduces to 


amy: _ ned 
pea =nkT = “a, (2-8) 





where D and M are the diffusion and mobility coefficients respectively. 
The last tenn is derived from the preceeding ones by the use of 


Einstein's relations 
~~ o = Se: (2+9) 


This generalization (equation (2-8)) is valid only for a plasma which 
is characterized by a small drift velocity, 

il <l¥. ( 2-10) 
The lest term in equation (2-3) is frequently written 


ol 


-Whv = - a Cam 
which defines the mobility, .- In the case of ioneneutral particle 
interactions, Y is a constant because the forces are inverse fifth 
power polarization forces. In all other cases we must introduce an 
appropriate average value, ¥ , which we will consicer independent 
of Vee 

One important consequence of the above assumptions is that the 
ion and electron mobilities and diffusion coefficients are scalar 
constants. 

Substitution of equations (2-8) and (2-11) into the momentum 


equation (2-3), and rearranging terms, yields, for the positive 


ions of mass Ms 


af oe sd ——. aS - — a 
Male + (Din) - mpeE —pl) xB = -[f, (Cae) 
ané for the electrons (mass m)>; 
mu-3t. 4 Geom. z foe. 
nM + VPM) + nk + wlixB = - FR, (2-13) 


in 





For stationary (steadyestate) situations in the plasma, we find 


that these equations reduce to 


FF yelexO = Fpune-E -Ds Vos, (2=2h) 


B. Current Density, Field and Potential Equations 
From Poisson's equation we have 


VE = & = Gms (2-15) 


D 


and from Maxwell's equations for the stationary case: 

VXE =O, (2416) 

B me —_— 

ae - ) 

Vx e(T" i), (2-17) 
These three equations, with the continuity equation (2=6) and the 
simplified momentum conservation equation (2=1h) express the 
governing conditions for a plasma thich complies with the assumptions 
ennunerated to date. These will be solved for current equations and 
for general field and potential relations. 


From equation (2-6) it is readily deduced that 


[hell 4+(VxF) (2-18) 
in@hich (V x E ) is a vector which we shall now determine 
Let us Gefine three orthogonal directions in the plasma, one parallel 
to the magnetic field, designated by the subscript "||", one 
perpendicular to the field and to any surfaces parallel to the field, 


"\", and the third transverse to the field and to the perpendicular 


> wp 
~<a 

e 

a 





direction "T". In cylindrical coordinates, with a field parallel to 
the z-axis and cylindrical boundaries parallel to the field these 
would be the 23 r5 and Y— directions respectively. The ion and 
electron particle currents at the sheath edge are equal at all points 
for a steady-state plasma because they must arrive in equal numbers 
to recombine. For a geometry in which the insulator surfaces coicide 
with coordinate surfaces (the only geometries for vhich we will be 


able to obtain exact solutions for the particle density distributions), 


=P =P, (2-192) 


eg li, ( 2-19b) 
+4 cll 
at the sheath edges. We now impose the condition that 


Vin, = Yn-=E, = 0, ( 2-20) 
This requirement is met either by a plasma which has a rotational 
axis of symmetry parallel to the magnetic field or one in which the 
field lines coincide with one set of orthogonal coordinate lines anc 
in which the resulting distribution is a function of only one coordinate. 
The latter condition would be met in the vicinity of a plane surfece 
parallel to the magnetic field lines in an "infinite" plasma. Here, 
for a y=z plane and a field parallel to the zeaxis, the perpendicular 
direction would be the x-direction and the transverse direction the 
y-direction. If these symmetry conditions are imposed, it is readily 
deduced that equations (2-19) apply throughout the plasma if the 
walls are insulators. If, however, the end surfaces are concuctors 
(i.e., electrodes), equation (2-19b) does not apply. The paraliel 

6 





current is then a conduction current, not a diffusion current, and is 
not treated by the developments to follow. This current does not 
affect the perpendicular currents, however, and such a plasma vill 


be treated as an infinite plasma in the longitudinal direction, with 


arr = Oo: (221) 
From equation (2=1h) and the symmetry conditions (2=20), it is 


immediately seen that the transverse component equation is: 


er Seal ee (2-22) 
From this equation, and equations (2-19), we can see that (VxF) 


of equation (2-18) has only a transverse canponent: 


Tee) i/o M : 
(VxF) =(VxF), = -(lt SE) i 4 (2023) 
The parallel and perpendicular component equations derived from 
equation (2=1l) are 
ls) - a ae Ee a DsiVinz (2-2) 


anc 


a = pyle = TUM EY SD (2-25) 


A simpie solution is possible if the assumption ef proportionality 


holds: 


War x Vn- = Vin ( 2226) 
Ne n- n° 


By eliminating E,, from equation (2-2), we obtain 


Sim eet Ne sind e=2 
F a DsViin- = - Dyn, Yan- es 
vhere Dray is a quasi-ambipolar diffusion coefficient for diffusion 


7 








parallel to a magnetic field and Ds is the ambipolar diffusion 
coefficient which is commonly developed for a piasma when B = O by 


infering congruence, 


ca en consent 


[’ a ( 2=26) 


from the current density divergence equality equation (2-6) and from 
the proportion: lity condition (2=26) (27 3. Dg is related to 


Schotthy's ambipolar diffusion coefficient, 


D+ pA eet DM 


J ( 2-29) 
Da. M, + MM. 


by the equations 


ioe (rt f-\Ny 
” ayn, + Men.) 


D.C is as), ( 2=30a) 


In equation (2-30a) we have employed the following relations for 


De ( 2-30) 


il 


l 


charge density and conductivity 


fp - os —7-)2 ( 2=31a) 
and 
\ 
Sean Pe, aT 
By substituting equation (2-22) into equation (2=25), we obtain 


Mey at ieee. = tMe Nz E, SD eae. ( 2~32) 


By the use of equation (2-19a), E, may be eliminated from the 


electron equation, to yield 


Mtn +4 le Ne + le WeN- t f- “ls y.) B ~ Dip t Dy, V ov (233) 
Me P-Ne N- i fle 





By simplifying and solving for ri » we obtain 


Zs Dwi 
h, = (Ut _ ty) Z 
i+) + CM ay B* 


wnere Dm) is a quasSieambipolar diffusion coefficient for diffusion 


= -Dm Vin, ( 2=31:) 


perpendicular to a magnetic field. Thus diffusion perpendicular 
to the magnetic field is reduced by an inverse dependence on the 
square of the magnetic field intensity. 
If there is no superimposed electric field, so that E = Es, 
where Es is the space-charge field, Es, may be obtained by substitt cing 
equation (2=3l)) into (2=32): 


Es, A. er (2—35) 
Similarly, 
FE, - Damo B- Van . Ds=D- Win (2236) 
NN fA— vy) M- Y) 


The transverse space charge field is zero aS a consequence of our 
assumpt ions that there is no transverse density gradient. 
For a quaSieneutral plasma, in which 
WANN, (2-37) 
it follows from equation (2=30) that Ds = Da, and that 


Fs, ——> Se un (2-38) 


N,N fy t p- Y) 
In the absence of a ragnetic field, equation (2-35) when sinpiified 
and added to equation (2-36), yields the single equation for the 
space=charge field 


= Ds- D. Vr ~ V 
ee at eG _ ———> Deer Vn R% 
/- oa M+tM- 1 * (eo? 











From the field equations, we obtain the following equation for 
‘ 
the space charge potential for a plasma with only a perpendicular 


density gradient: 


D--(I+ B Dm, S 
Vi = = tg a, (21:0) 


In the absence of a magnetic field, equation (2-39) yields 


oa - 2 
V = Aa | Gee ee) OY 


The zero for potential in both of these equations is taken as the 


"center" of the plasmas that point (or region) where the density 


gradient is zero. 


Cs fhe Diffusion Tensor 
The current density equations (2-22), (2-27) and (23h) can be 


combined into one equation 


—_ ao 
& = —D im. Vn, (2212) 


which, in rectangular coordinates, is 


: Dm, Dre ven 
eo = Deas 4 Dm =) Vay ¥). (2=1;3) 


[+s I O O Dru, Vi NE 


) 


and in cylindrical coordinates is 


ee Dn, Dae O Vn 
[he One Dray O O ( 2=))1) 


[tz O) © Dm, | Von ts 
10 


\\ 








= 
- 


The components of the ciffusion tensor are 


Den, i. Ds, (2=lta) 

Dm, = ane OM (2=)£b) 
t+[t + Pn L) Aye BY ) 

Dey = f+ EB Dey (2=l;5c) 

Bs - 8h (2-82) 


We shail now examine each of these components. Equation 
(2-):5a) reflects the fact that the magnetic field does not affect 
diffusion in the direction parallel to itself. Further, we have 
shown that 
Dra, =). pao > Vo. é (2=1:6a) 
pngenyigeewt 
Allis [2] has also shown that 


Do Se 


On = CO 


(2=L6b) 


) 


where g; is the conductivity at the center of the plasma. The quasi-= 
neutral limit (3-lj}6a) is applicable to a plasma in which the 
mobilities (and consequently the conductivity) are so great that 
negligible space charge can develop. The above limits are the Low 
and high limits for the ambipolar diffusion coefficient. 

The perpendicular component is more complex. First, as the 
magnetic field approaches zero, 


Dr, ——> Ds (2-h7a) 


t+ 68-0 


11 








As B becomes very large, Dm) approaches zero as pe 


mere em eee) ee ( 22):7b} 
oa Boece [4s p- Bs B—>co . 


From equation (2-l6a) it is seen that 


—- _— fy (2-h7c) 
“pze Tap, 


and from (2-60a) that 
J euie 
es —— fe ( 2-h7d) 
Ono (d(KE- Ar) pf B 

which can be non-zero only if (oo O. As shown in equation (2=55), 
for a long, narrow plasma this can be true only if 

bre a 
L G70 / — fl B~ . 


From these equations it is apparent that for the parallel and 


(2:7e) 


perpendicular diffusion tensor components in a magnetic field the high 
and low limits are the corresponding electron diffusion values or the 
quaSi-neutral ambipolar values respectively. 


For the transverse components, we obtain 





Dm.r — ie : . ( 2=1:8a) 
/ -f- Ly ae ae ee (3 

EOE ( 2=1;8b) 

ly 2 H- p+ ime 

aon 
2 Se > 215 ce ena 

B70 " hts 15 Boe @, (2-Ltc) 
B eG =M2Dsb Bo - Base 


le 








The maximum transverse internal current will flow when Dm.+- area 
maximum. For the quasi-neutral plasma, equation (2-l8b), this 


occurs when the magnetic field has a magnitude given by 


\ 


B= Wiest ( 2=1:9 
For helium, using data from reference [ds » this is approximately 0.12 
weber /n@, The magnitude of the current density resulting from this 
field can be determined from equation (2=22) 


- 


Yn 


C(jA tp) DY LA 
= [- re ra (' a — — ee (> 2_50 
ax e( ey a wT ad tft , a 9 ae ( 2 ) 


where . is a unit vector in the transverse direction. The magnetic 
field generated by this current will oppose the applied field, and 
accounts for the observed diamagnetism of the plasma. This is the 
term which satisfies equation (2-17), and it must be small compared to 
B to justify neglecting its effects on the other plasma equations. 

In the quasi-neutral plasma, with 71, 2?i_, the diffusion tensor 
reduces to 
+ MDB 
| +p1,fa-3* 


Da O (2=51) 
| + [As p- B* 


O es 





D. The Continuity Equation 
By substituting equation (2-2) into the continuity equation (2-6), 


ve obtain 


en ~V:(De, Vn) = YY. (252) 





a second order partial differential equation. For the steady state, 


this reduces to 


V(Dn Vn.) = Yn (2-53) 
whereas for a decaying plasma (for which ) = 0), it becomes the | 


diffusion equation 


vs = V-(DmsVn é.) =O (2~5h:) 


In order to simplify equation (2-5), we must examine its com 
ponents. The components of the diffusion tensor, equations (2=l5), 
are constants because of our assumptions and consequently the 
divergence of the diffusion tensor is Zero. In the absence of a 
magnetic fields D,,is the scalar D,, and equation (2=52) becomes the 
familiar continuity equations: 


on 


Ds V7n. + Vn = SE (2-55) 


For the steady state, we can further simplify this by employing the 


assumption of proportionality, equation (2-26), to obtain 


. 
oo =~ (2-56) 


In the presence of a magnetic field, we shail limit ourselves to 





orthogonal geometries in which equation (2-53) can be separated into 


[Ora Vin. sie Dm, Vine |- PH. =O (2-57) 
In equation (2=57) we have made explicit use of the requirement that 
the transverse density gradient be zero. By employing the assunp= 
tion of proportionality (equation (2=26)) and a theorem from vector 


analysis, we see that 


Ly 





ie 4 
with a similar oe for the perpendicular components. Thus we 


Y n- mt SSS. (_ wu 7 (Yi) (2-48) 


can divide equation (2-57) by n_ and expand it to obtain 


9 
7 


3 .( Vr V V, Via)” 
ON (2) + Digi) D(H) + Dy B) con 


If we then define a vector u such that 


tA = Upt U,, = ao da aoe = yy (2= 
we obtain 


7) oo —— aot 2 2 
Dra, ea Dra, Vn eee Dra Uy & Dm Wyo4Y=O, (2-61) 
which is a first order partial differential equation in u. In the 
absence of a magnetic field, D is a scalar (D.); and equation (2-55) 


can be sinilarly expressed (in any coordinate system) as 


— 
ae erD 


J > i VY 
Vel + +B =O. =) 


The significance of the transformation of this equation from a 
second order differential equation in particle density, n, toa 
first order differential equation in u will become apparent in the 


discussion of boundary conditions in a later section. 


15 





E. Boundary Conditions 
1. Derivation 

We shall now develop an equation for the boundary condition 
and evaluate it for a few simplified cases. We consider a quasi~ 
neutral stationary plasma in contact with a plane wall. The mag- 
netic field, if any, is parallel to the wall. 

There are two conditions which we must fulfill. The first 
is that the electron and ion currents to the insulator surface 
must be equals the other is that the current must be continuous at 
the limit of the diffusion region, which can be at the edge of the 
sheath or at the edge of the free-fall Zone. 

We have to distinguish different sheath theories accord= 


ing to whether D_, > eee DS or D-1*D,, , where 
D = (2-66) 


a.) The case Day >> Diag, 16m current anertia—limeecd. 





Let us first examine the problem in which D_, >> Di, 
In this case a negative wall charge builds up until the electron 
current is reduced to equality with the positive ion current. 


The electrons will be treated as having a Maxwellian distribution 


ee (2-67) 
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where the potential V at the sheath edge is taken as Zero, -e is the 
charge of the electron, and the subscript "o" means evaluation at the 
sheath edge. This notation will be used throughout this section. 

Ve find that two possibilities exist for the ion motion through 
the sheath: (1) their motion may be inertia=limited, or (2) their 
motion may be mobility=-limited. In the inertia-limited case, the 
sheath is a free-fall region and the current in the sheath is that 


given by kinetic theory: 


inert 7 Vi : 
to a ? (2-68) 
We must next evaluate the diffusion current to the sheath 
edge, which is 
a 
i Ror =", 
s Vin, . (2-59) 


oe [+ A4p-B~ 
If the electrons and ions are in equilibrium within the body of the 


plasma 


a CUP ay | (2~69a) 
In order that the number of particles varried to the sheath by 
the above diffusion current equal the number carried to the wall as 
expressed by equation (2-68), the normal component of the diffusion 


current must equal the sheath current. Therefore, we have the 





equation 
a me poe 
Mev 3), v7.) Vine (2-70) 
4 / ie : 


From this equation we see that 





Yn of 
(4 O - Me a A+, a7!) 


a result obtained by Ecker [3], where, for this case 


ref 





ro: -B*). (2~70a) 
It will be seen that the boundary condition can always be ex- 
pressed in the form (2-71), with the value of & varying with the 
sheath characteristics. Let us investigate the potential function V 
for this condition. If there is no magnetic field, we see from the 


conservation of energy that 


2 ara 
Cine cl? ee (2-72) 
z Zr 


where M is the ionic mass. From this we can derive the following 


equation for the ion current density: 


/ V (2-73) 
Ay = UCM, ~mCV,, = CV. / —F> 


Let us now employ a coordinate system in which the x-direction 
is normal to the sheath and positive towards the wall, and whose 
origin is at the sheath edge. Solving equations (2-67) and (2-73) 
for n. and ny respectively, and substituting in Poisson's equation 


(2-15), we find 


2 V | 

AV eC C7 C 

— ae = ae - aoe herd <. ane 

dx* = (n- n,) E, € ae = (22710 
Multiplying by AV. and integrating yields: RT, 

Ax 

Av \2 2n,\eT. eV oF { > 
S| = = 41 ee als eee = 
a CE. (. Rt) + =| sn ) + Eo, (2-75) 
wnere we have set | = =f, and Vv} = QO. This result has been 


reported by Ecker{5]. For the sheath to be stable, we -E must 


increase with x. Therefore the quantity in the brackets of equation 
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(2-75) must be positive. Expanding this quantity in series, we find 


ai V 
eb oer + <i (- Zur, 7 reek = (2~76) 
from which we see that 
| 
eae (2-77) 


is the condition for the sheath to be stable. Bohm [7] has calculated 
the function V for two values of T, by numerical integration. 
If we neglect the electron density in Poisson's equation, and 


integrate, we obtain the Langmuir relation (8): 


‘Ls 
V = -(2i)5 (44x) (2-78) 
pee tne case =). >> Der sien current iN ~limited 
We shall next examine the case in which the ion current 
through the sheath is mobility limited; i.e., the ions move in the 
field of the necative wall charge according to their mobility. Now 
the electron current density striking the vall is obtained from kinetic 


theory and the electron density relation (2-67): 
, eV.) 


Las = "ere 
ee ON : 
‘ a . ut (2-79) 


where V(d_) is evaluated at the edge of the electron free=flight zone 
(suitably corrected to account for curvature of the electron paths in 
the magnetic field). Again equating this current with the diffusion 


current brought to the sheath, we obtain: 


e Vd-) 
a ., 2 
jy. om 20.) Van 
a = = (2-80) 
tT | ae Bo. 


Ie 





In order to further simplify this relation for the case in which the 


electrons and ions are in equilibrium, we use the relation 


= 2 — . 
mv, = Mv, (2-81) 


from which we See that 


= eV(d-) 
oe - #7, we s (| + afl B). ae) 


Thus we obtain for the parameter cé : 


eV.) 

= 7 e ATE ‘Cea BD. (2-83) 
In order to evaluate this expression, it 1s necessary to know the 
potential function V(x). The ion current density in the sheath 
(neglecting the magnetic field) is: 

7 AV 

j+ — cea f4, V1,° ie = — fine Ax ; (2=8)1) 
where 
( ( 2=-8ha) 


5 tue aos 
jie pepe B* 


Again employing Poisson's SCUaiOn and solving equations (2-8):) and 


(2-67) for n+ and no respectively, we obtain: 





> V 
dV C\ ye. , —t (2-852) 


Ax* ~— & pine dt, | 
| V : 
$e) et y be (2-85) 
CAL e 


V 
tS = QO, 


GEw); 


At the sheath edge (where x = 0, V = 0, and a 


a result consistent with the assumption of quasi-neutrality at the 
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sheath boundary. Equation (2=-85b) can be integrated once, and yields 


a first order differential equation for Vs: 


(GE) = Spefeex + er (e% -)) + Bt ai 


AIX 


Further integration, to obtain V(x), must be performed numerically. 
However, if we again neglect the electran density, equation (2-86) 


can be integrated to yields: 


3 4,x° 
Ve [3 (2-87) 


e. Ihe case D.,. ~ De: 
An approximation for the case D., ~ D+, can be easily 
obtained for inertia-limited ion current. By neglecting 4, B* 


as compared to 1 in the equality 


ith sree yin ee 

[eae | +p B™ , Meare) 
we obtain 

= zx J+ 2B, (2-89) 


meom vnaen it iseseen that 


Pe ee 4 
M- = Mf B (2-90) 


Employing this relation in equation (2-70a) yields 


Ke (+E = us) (2-91) 


If we neglect the last term in the parenthesis as compared with one, 


al 





and employ Einstein's relation 


D, kh 





a. (2-92) 
Mb. Cr 
we obtain 
yee T 
a= g(1+ =), (2~93) 


ad. The case D4pD—) 


For the plasma in which D4, D., , a positive wall charge 
must build up until the ion and electron currents are equal. The ions, 


assuming Maxvellian distribution, will have a density within the sheath 


given by 
Th, = Mo€ ) (2=9)) 


where the ionic charge is +e. 


The ion current to the wall is obtained from kinetic theory as 


_ eVia.) 
NN 4 ee 
L} d. - Lf ) 


where the evaluation of V(dy) is at an ionic mean-free path, d,, 





(227) 


which has been adjusted to include the effects of the magnetic field. 
Equating this to the diffusion current to the sheath yields the 
boundary condition (2-71) in which 
. 3 _ eV(ds) : 
aera BT (+ 5-8) (2-96) 
In order to evaluate the function V, it is necessary to employ 
the equation for the electron current (neglecting the effects of 
the magnetic field): 


aie 





ey 


‘ AV 20 
qd” = — fle MoE E.. = = fly ree = fAl-ii- Ee ax Ps ce 97) 
By the same procedure used for the mobility=limited ions when 


D-,>? Ds) 5 we obtain 


(vy _ eel iene Ti(o" -)k i ae 
dx Ey Z e ; (2-98) 
2. Evaluation of & 

The overeall prohlem of evaluating “ of the boundary condi- 
tion equation (2-71) is a much more complicated problem than indica- 
ted hy the derivations of the previous section. The influence of the 
transition region, the effect of non-planar geometries, deviation 
from thermal equilibrium within the plasma proper (and in the trensi- 
tion region), and secondary emission are just a few of the factors 
which can affect this quantity. The transition region, which joins 
in a continuous fashion the sheath and the main body of the plasma, 
presents particularly difficult problems which are beyond the scope 
of this thesis. In the computation of & for a surface whose dimen- 
sicns are comparable to or less than the mean free path, it is 
especially important to compute the fraction of particles brought 
to the sheath edge which would miss the surface in their flight 
through the sheath and would pass back into the plasma. Thus, for 
the limit where the surface is very much snaller than the mean free 
path, as for a Langmuir probe, & can become vanishingly small. 

3. Discussion 
An understanding of the processes occuring vithin the sheath 


has been shown to be necessary for the evaluation of the boundary 


ae 





condition equation (2-71). We have seen that the requirement that 

the current be continuous results in this condition which must be 

met at the surfece of the ion free-fall zone (the region within one 

ionic mean free path, ), , of the wall) or at the surface of the 

sheath, whichever is farther from the vall. As the sheath thickness 

and A, become very much smaller than other dimensions, equation | 


(2-71) may be replaced by its limit for vanishing A, : | 


11) att = 0. (ae) | 


This is the generally=used houndary condition, but for small values 


of « and large ),, the error in (2-99) can be appreciable. From 


(2=71) we obtain: 


nN) on % Mt mMn = ny (1-0) (2=100) 


where the subscript d, refers to evaluation at a distance A, from the 
wall. 
The importance of equation (2=71) and of the transformation of 

the continuity equations to the forms (2-61) and (2-62) is now 

readily seen. The normal component of the vector u defined by 
equation (2-60) must have a constant magnitude, =e along the sheath 
edge. If there e more than one wall, anc the walls involved are 
orthogonal, the density gradients establishec to satisfy this condition 
al each sheath edce will not interfere with each other. Parallel 
surfaces, on the other hand, create an eigenvalue problem beceuse u 


must satisfy the boundary condition at one sheath edge, must change 
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as prescribed by the continuity equation (2-61) as we move across 

the plasma, and musi acain satisfy uniquely the bouncary condition 

at the other surface. If the walls are neither orthogonal nor parallel, 
the density gradients established to meet the boundary condition for 

one surface will have along the other surface normal components which 
are not constant, and vice versa. This will result in a violation 
of equation (2=71). Consequently, for non=-orthogonal geometries we 
are unable to predict the steady state distributions with the equations 
developed in this thesis except by an approximate method to be 
discussed in a later section. 
eecat Losses 

The only source of heat loss considered in this thesis is that which 

occurs at an insulated surface as a consequence of the particles 
impinging on the surfece and giving up their enercy. This energy can 
come from several sources: (1) the recombination energy, eV;, fron 
each electron-ion recombination, (2) the kinetic energy of the particles, 
and (3) excitation energy of excited ions which return to their stable 
state upon colliding with a wall. All of these are proportional 

to the particle current density to the surfaces. As a simplification 
for illustrative purposes in later calculations we shall limit ourselves 
to plasmas whose temperatures are sufficiently low (below about 10° 
degrees) that only the recombination and electron kinetic energies 


are appreciable. The heat Liberated at a surface is then 


As 
Eun = Gf (evi+ S), (2-101) 


where the subscript "n" means normal to the surface and foi) aS 


expressed in electron volts. 
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SOLUTIONS OF THE CONTINUITY EQUATION 


A. General Solutions 

This section discusses a procedure to obtain particle distri- 
bution functions in several coordinate systems, and presents some 
of the general features of the solutions. The following sections 
develop many solutions for specific geometries. 

1. Exact Solutions 

ae General Developments 
With no magnetic field, we have exact solutions in the 

eleven orthogonal coordinate systems in which the wave equation is 
separable when the boundary surfaces coincide with coordinate planes. 
The procedure we shall describe applies to the non-steady state 


plasma for which we repeat the general continuity equation (2-55) 


2 
OV. +n. = 3b, — 


wnich becomes a wave equation for the steady-state 


- yy 
VO em 3 j (aage) 


for a decaying plasma with no ionization it is 


2 on 
D; V mM.= >t (3=1b) 


with a steady-state counterpart (Laplace's equation) 


Vin. = O. (3~ic) 


In the presence of a magnetic field, we further require that the 
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coordinate system be one in which the continuity equation may be ex-= 


panded as 


Dm, Vi 2. + Dm, Vin. + vy-< an 


mm (322) 


(see equation (2-57), 

Equations (3-1) and (3-2) are separatle in the geometrical coord= 
inates and time, and the general solution for a generalized coordinate 
Systems U, Vy Wy, may be expressed as 

n(u,v,w,t) = N(u,v,w,)T(t), ( 3~3) 
Substitution of equation (3=3) in equation (3-2), and division by NT 
yields 
2 
Dm Vi, MV a Dm, Vi MV 7 Ta L 


NV 7. V = ia = t, (32) 


where =k, iS a Separation constant. The time solution is 





i a G on eit (3=5) 
and then 
Dyn, VN Dm, VN (3-6) 
ay) + J = = (le, + y) 


In the absence of a magnetic field, equation (3-6) reduces to 


Z 
A/ 
ay al le”, (3+7) 


in which 
b +v 
z y 7a 
L. = ‘oe ; (3 7a) 
and, for the steady state (ky = 0) 
a ae : 
bbe = be. (3=7b) 
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Our approach to the sclution of equations (3=6) and (2-7) is to 
employ the technique of separation of variables. Thus we assume 


that the solution is of the form 


N(u,vyw) = £ (u)fo(v) £3), (3-8) 

and, by substitution into equation (3-6) or (3=7) (as appropriate), 
we obtain three differential equations, one in terms of each of the 

coordinate variables, which the functions f;, fo, and fg3must separately 
satisfy. Magnus and Oberhettinger [9] have summarized the results for 
eight coordinate systems. Because the separated equations are of 
the second order, there will be two linearly independent sets of 
eigenfunctions for each coordinate, FL) (x, yu) and F(x, eu) 5 
GD kn ov) and G2) (ke ov), and HO) (ke gw) and Hh) (kaw), which satisfy 
the differential equation in the corresponding variable. We shall cal} 
the functions with the superscript (1) functions of the first Kind 
and the others functions of the second kind. The most general solue 
tion to these equations will then be a sum of the products of all 


Such linearly independent solutions of the form 


Nona lt,¥, 0) = Cag F “egu) + CogF% be, LC .C"Cep v) HCme( a lX 
a L Con H” (eg, w)+ Can H ew) | - (3=9) 


The complete soiution for the particle density distribution is 


then 
Play, vit) = FE Norn (“rd TG) ae 
-Kt 
= d Ai,,,(uya)o™ (3-10a) 
4m, : a 
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The summation is over all possible values of [, m, and n. Ky, kins 
and k, are eigenvalues or spearation constants (depending on the 
coordinate system chosen) corresponding to f, m, and n respectively. 
We can sum over only those eigenfunctions which belong to a degenerated 
eigenvalue (as determined by the boundary conditions). It should be 
noted that the time=depencent solutions to be developed herein are 
for those problems for which the boundary condition (2-71) or (2-99) 
ig still applicable. A plasma decaying from a steady-state distri- 
bution upon cessation of ionization is considered to be such a 
plasma. 
b. Cartesian Coordinates 
An example of the above method of solution is the solution 

in certesian coordinates in which the functions of the first and 
second kinds may all be cosines and sines respectively, or in which 


one or two pairs of the functions may be exponentials. For example: 


Ni, 4 2) = ‘mg Lay coal t,x) Cx rin.(k ALG, cork, 4) + Goa Ny 


( 3=11) 
XC pcealla?) + Cy, orelkye)|, 
By substitution of this selution in equation (3-7) , we see that 
Kk? of equations (3=-7a) and (3-7b) is given by 
2 eee > a 
2 2. (b+ k,2+k,*), (3+112) 
yn 


and consequentiy that the time solution (equation (3=5)) is 


cekt cg tbh cum 


Ia, & - Cie = 


2g 








Obviously, when 

LL? = (b? b 2 b*) - _ (3-12) 

SUES tKy thy )= 5, 
Pa 
the time dependence vanishes, and the plasma is "stationary." 
Another possible solution, involving exponentials, is 
box ay | by aC { | 

NOs 2) = 2 tue rly, 6 WG, € 4G e ] G, cole,2) Hy, s (U2), (3-13) 


where Ke is given by ke = 2 adele st, ), (3=-13a) 
ym 


c. Cylindrical Coordirates 
There is an independent eigenvalue, k;, for each set of 


functions only in cartesian ccordinates. For any other coordinate 


system u,v,w, defined by a transformation equation of the form 


x = x(u,V,w) 
y = y(u,v,w) (3-1h) 
Zs Z2(U;V5w)) 


those eigenfunctions of the variables u,v,w which appear together in 
any of equations (3-1) will be connected through a separation variable 
which is required to reduce equation (3-h) or (3=7) to separate 
equations in each of the coordinates, Each such separation variable 
required for the above purpose reduces the number of eigenvalues 
available for the satisfucticn of the boundary conditions. For exs= 


ample, for cylindrical coordinates, where 


x = rcoll 
y =r ew 
Z22 (3-15) 


a 


the eigenfuncticns in r and (? are related by the separation constant, 


m, aS indicated in the weli-slnown solution 
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M(r, 4,2) my <. Rér) E(®)Z(z) (3-16) 


where 
Riv) a Cae, Sab kr) + Caen Nin Cer) ke, 40 
= Cc + a. k,=0 
= Cher + Cs <0 54ee) 
SC) ee i A ; 
Z(t) =CeyCoalkzz) +Cy, sam_(ty2) wn : 
and = Cz a C7 n ») 
ihn - (4 bed (3-16b) 


An alternative Z(z) solution is 


b., sens 
SAEs) = Cr € “4 Cae, ke, 40 (3=16c) 
with 
Be 5b bn (3-162) 


Jn(k r) and N,(k r) are Bessel functions of order m of the first and 


second kinds respectively (Bessel and Neumann furctions). 


Cc. Spherical, Prolate anc Oblate Ellipsoical Coordinates 


In the other coordinate systems to be discussed, the general 


solutions of (3-7) are: 


Spherical: 


M59?) = 2 R(r)\@C) £(¢), ward 
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where 


Rir) = com (elt r) FG, (tr) L #0 
See Tata E, = 0 


EX) = Cu, "7 aoe) +Cs4Gl wv (ca20) 


Sly) = Cm conta) + C5, Sun (re yf) mMt0 (3-172) 
= CY a C7 m=O, 
and 
pr. g§ va (3=17b) 


and in which Py( cose) and Gn(cos®) are Legendre Polynomials of the first 
and second kinds respectively and j,(k r) and n,(k r) are spherical 
Bessel functions of order n of the first and second kinds respectively. 


The latter are defined by the equations 
TT 
Aun (ber) = 2ky Tet (ee) 
Ayn (by r) = {= Nope (be 1) 


Prolate Ellipsoidal (See Appendix I for diagram and transformation 


(3=17c) 


equations: notation empicyed is that of Flammer 6) with slight modie- 


fication): 
(1) (2) 
Non 84) = Caran Sree Ca) + Catan Somn Cod] * 
a nae? (4,3) Cae Ro (6,8) Kenatoeet)t (4, Sint ee 


where 


(3-362) 


e 
RFC, 


and d is the interfocal distance. San) (6) 9) and sf2)(c, %) are spheroidal 
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angle functions of the first and second kinds and RAM (cy, =) and 
e,) g ) are spheroidal radial functions. See Appendix II for 


tables and graphs of these functions. 

Oblate Ellipsoidal functions may be obtained from the prolate func- 
tions by the substitution of wic, for Cy where appearing in equation 
(3-18). Equation (3=-18a) is unchanged by the transformation. 

Similar equations are obtained in parabolic cylinder coordinates 
(where the functions include Weber-sHermite functions), parabolic 
coordinates (Whittaker functions), and elliptic cylinder coordinates 
(Mathieu functions). These are not developed further in this thesis 
because their geometries are not useful in the problems with which we 
shall deal. 

d. Solutions in the Presence of a Magnetic Field 

In the presence of a magnetic field, Dm is a tensor, and 
equation (3-6) must be employed. In order to make a separation of 
the form (3-6), we are limited to cartesian coordinates or cylindrical 
coordinates for the problems which we shall treat. Parabolic cylinder 
and elliptic cylinder coordinates are also suitable, but we shall not 
need them for the geometries with which we shall be concerned. Firther, 
if we consider the assumptions under thich Dn was developed, we may 
immediately limit the treatment to the case where m and C are Zero, 
Since no gradients in the ¥ direction can exist. 

In the coordimate systems selected, equation (3-5) can 


be reduced by separation of variables, by assuming 


N =N,WN,, (3=19) 


Bo 





~- ——$—— ree 


and thus obtaining two equations 


WN, a = ie. be 


co, =. a 2 -k, (3=20a) 
ai J 
and Ma 
z o 2. 
nee = 2 = —k, (3-20) 
a Oia 
where 
Dk, + Dw ky = R? = kV, (3=20c) 


The solution to the above equations is the same as in the absence of a 
field, except that equations (3-12) and (3-16d) become (for a field in 


the zedirection): 


ee b+ Y= E Dre, (ey kw) * Den, ea" : (3-12") 


Liem > 


and 


| , 2 (3-15d? ) 
kt =e kV Z po i Dra, Fn \ : 


ee Further Simplification of Previous Solutions 
The equations developed so far can be simplified considerably 
by considering some of the restrictions upon the solutions. The require= 
ment of a single=-valuedness linits J, m, and nto be integers, and it 
is possible to further limit them to positive values without loss of 
generality. Further, if the plasma possesses rotational symmetry 
(which it must in the presence of a magnetic field under our assump- 


tions), then n and C6 are both zeco in the equations for N(u,v,w) 
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(with the exception of rectangular coordinates). The requirement 
that the density be non=negative everywhere within a given region will 
further limit the solutions. Because all the functions discussed 
(except the exponential function) are periodic, this latter requirement 
usually limits the steady-state solutions for bounded regions to the 
first positive segment of the applicable function, and values of k,, 
kn» and/or ky, are thus determined by the dimensions of the plasma 
container. From the boundary condition (2-71), the physical reason for 
this effect is apparent. u oY has a definite value at the edge of 
the sheath. paGabion (3-7) can be transforned into an equation in u 


in the same manner as in the derivation of equation (2=62), yielding 
Ve + me = he (3=21) 


It is the value of k* in this equation that determines the "trajece= 
Eory" of the vector U across the plasma to the sheath cn the "opposite" 
side where u must again satisfy the boundary condition (2-71). Because 
equation (3-22) is a first order equation in u, the problem is completely 
stated (except for a condition to determine an arbitrary normalization 
constant), and only one value of ke is possible. 

Further, and obvious, simplification can be accomplished for many 
of the solutions stated herein if they must describe a plasma in a 
region in thich some of the eigenfunctions are singular. In parti- 


cular, the following functions have singularities for the values indicated: 


Nm£0), Mn(0), O21), SkXc,21), S42) (-ic, 41). 
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2. Approximate Solutions 


There are a limited number of geometries for which a stable 
ambipolar solution may be obtained by the separation procedure de- 
scribed herein. It is, however, possible to obtain an approximate 
solution for the particle density distribution if the geometry is 
not too radically different from a geometry for which an exact 
solution is dmtainable. 

In the discussion of the boundary conditions (Section II-E=3), 
it was shown that the conditions at the sheath edge control the den- 
sity distribution in the vicinity of a surface. If there are two or 
more surfaces bounding the plasma in Such a manner that it is not 
possible to obtain an exact solution for the particle density dis- 
tribution, then we divide the plasma into "Zones of influence" 
about each of these surfaces. We then calculate the density dis- 
tribution for each zone as determined by the boundary conditions at 
the insulator surfaces belonging to that zone. In the region near 
such surfaces we expect our solution to be very nearly correct. There 
will be a transition region between the zones of influence in which 
we expect the greatest difference between the model and the actual 
plasma. In this region there is either a departure from the ambi- 
polar diffusion or a time dependence (non-stationary plasma), or 
both. However, we expect that processes in this transition 
region do not affect the total solutions significantly in many 


problems. A measure of the importance 
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of these processes is the degree to which the solutions for the 
different zones conflict in the transition region. When the dis-~ 
continuity across the transition region is small, we expect that our 
solution correctly describes the plasma in the immediate vicinity 

of the insulator walls and describes it approximately (within an 
order of magnitude) even in the transition region. 

In order to discuss this approximation more quantitatively, let 
us consider two walls, W; and Wy, bounding a partly-ionized gas, and 
assume that the configurations are such that the model described in 
this thesis can be applied to each wall in the absence of the other. 
Thus, two distribution functions, n, (q;) and no(q;), are obtained 
(with q; being the coordinates). If, through adjustment of the 
spearation variables of equations (3-6) and (3-7) and of the normali- 
zation constants, it is possible to make 


Ni, = “2 (3-2ha) 


and 


Er , oe (3=2))b) 


everywhere, then the solution is the exact solution. If this cannot 
be done, then we must find a surface, S, which separates the zones 
of influence of the walls W, and Woe The separate solutions for 
these zones, joined at S, form the zero-order approximation of our 
solution. 

There are three general procedures by which we might develop a 
criterion for determining the position of the surface S. We shall 
briefly discuss all three in the following paragraphs, but shall not 
attempt a rigorous proof. 
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One procedure is to employ a perturbation theory in which we 
define 
| 
M2 Not We + +. (9=25) 


and 


F = E. o Ye ae (3-26) 


where vi is a parameter defined by the equation 


ie Vee op (Mon-) (3-27) 


in which n,, is defined as the maximum density in the plasma. From 


equation (3-27) and Poisson's equation it is apparent that 


as fo » 
M a (3-28) 


For our zero approximation, we impose the conditions: 
Vilvewnene (3-29) 


and 


ee 7 ie - Le (3=30) 
in each zone of influence. 
In order to determine a criterion for positioning the surface 
between the zones of influence, S, (i.e., in order to find an adapted 
zero-order approximation), we must look at the first order terms. 


Neglecting quadratic terms in y » we obtain 


K a a (E. n,' + To aoe D, Va! (3-31a) 


and 


—- 


pe ~ —M(Eon- rné')- D.V7- (3=31b) 
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If we further assume that either there is no time dependence or that 
the time dependence can be accounted for by replacing all quantities 
in the following equations by their time-averages, we obtain the 


following from the continuity equation: 


— 


Vel +YN, t Walt aki st oon zy (aen(Fu, aa ) ( 3=32a) 
and 


Vs rt Bein! y 5 (Fer) A (ft (ies (oat + S6-F (Ti. je (3-32b) 


where the subscript "S' means evaluation at the surface S, $(7F-f.) 
is the Dirac delta function, and n is the unit normal to the surface 
Se The last two terms in the above equations represent particles 
which would have to be created or absorbed at the imaginary surface 
S in order to supply the currents indicated. The first order density 
terms can be either positive or negative. It is apparent from (3-32) 


that at the surface S, 
a = i — = 
ye pe P45.) - -(f,. - ee ee: = ( 3=33a) 


If the terms of equation (3-33) are identically zero, then the 
first order solution can be the'trivial solution 
Ns =O (3-3) 
Any deviation of the first order solution from equation (3-3) is 
obviously a function of the magnitude of the first term in equation 
(3=33). The surface S is thus the surface at which this term iS a 


minimum, If there were a surface along which 


(Dm Yn,,) A = (Dm Vor) y (3=35) 


oy 





everywhere, subject to the condition that the density solutions were 
also equal at the surface, then the op (Poi.- Me.) terms would 
vanish. This is not possible generally, however, and a least-mean- 
square procedure should be applied. Because this is not feasible 

for many problems, the following criterion for determining the 
position of S, applicable at a given point in which we are interested, 


is employed to obtain our adjusted Zero-order approximation: 


| Div Ving, cer, | = [Um Wo, cot, | (3-36) 


where Y, is the angle between Vn., and the surface normal and Yu 


between V/”,,. and the surface normal. Further, if 


Mn = 617) (3-37a) 
and 
Now * Cate (4 :) (3=37b) 
then at the surface S 
aie =VYboz5 * Ci, (4i)< = C.4,4.). (3-38) 


Solving for Cy, Substituting the result in equation (3~=36), and 


dividing both sides by C,f}(q;). yields 


oY a VEG 
Bm AS cat] = [Bane or] Go 
i q) S _ Js 


Because only order of magnitude results will be obtainable from 
this approximation, it is reasonable to assume as a first approxima- 
tion that 

coms, » cook, (3-0) 
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By employing equations (3-37), we obtain the following conditions 


defining the surface of separation 


fp Ve imams Vio | 
my Ney Wes Wor 


The second procedure to derive equation (3-1) is to consider 





(3~=)j1) 








the relative influence of each boundary at all points within the 
plasma, and define the surface S as that surface at which the in- 
fluences of the two boundaries are equal. Some measure of this 
influence must first be found. If we have in the neighborhood of 
a boundary a solution for the particle density which has a maximum, 
then certainly beyond that maximum the plasma is not influenced 
appreciably by that boundary. For the region between the surface 
and the maximum, the measure of a boundary's influence could be 
the current created in response to its presence. This procedure 
leads to equation (3-1) with an added restriction that if there 
are two points at which this equation is satisfied, we must select 
the one within the zones defined by the walls and the first maxima 
of their respective solutions whenever possible. 

The third procedure to determine the position of the surface 
S would be to calculate the probability that a particle at a certain 
position would go to the one wall or the other. The surface S 
dividing the two zones would be the surface along which particles 
would have equal probabilities of reaching all the walls in question. 
This might be calculated by stochastic methods or from the diffusion 


equations. We did not attempt the former. The latter method is 
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to find the surface from which the mathematical solution indicates 
that as many of the particles produced within a differential volume 
straddling the surface would go to the one zone as to the other. 


This procedure also leads to equation (3-l1) as the criterion. 
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Be Specific Solutions of the Continuity Equations. 


1. General Conditions 


In order to obtain an eStimate of the magnitude of certain 


quantities in the equations to follow, we shall use the following 


sample data, some of which is taken from Brown [lh] : 
gas: Helium 
p @= Imm Hg. 


T 700°K 


8 


I 
fi 


50,000°K 7.3 eV. 
ha *& 104 cm? /voitsec =o m? /volt=sec 
(ome 10° om* /voltesec = 75 n° /volt=sec 


a, me 06 n</s ec 


F 
R 





i, 
So 
} 


D kT 4, Ss x ~— 30 n*/sec 


De os D, = enh > Lely n@/sec 


~ hh volts 


x 
= 


Di- Ds 
fA tH - 
ev; = 24.58 eV 


et 
eV; + =e 29 eV 


E «= 300 volts/n 


S i Se VO eNO 


where "n' is the density in the center of the plasma or in an un= 


disturbed, constant=density region. 


(3-50) 


For the calculations of this section we assume that the mean free 


paths of the ions and electrons are both sufficiently smaller than any 


boundary dimensions involved, so that boundary condition (2-99), 


0% v* 


Vall ~ Oo, 
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is applicable. The more general boundary condition (2-71) will be 
employed in the discussion of cylindrical probes in the next chapter, 
and the method of employment of that boundary condition will be dis= 
cussed then. 
A quantity often quoted and employed in the literature is the 
diffusion length, A , which is defined by the equation 
Ps 
Rap ae Ee Ut. Ve. _ 
we shall compute this quantity for the various geometries in this 
sectione 
Finally, we shall calculate the "Zone of influence" of many of 
the surfaces to be discussed. This we define as that region in which 
the particle density is less than 50% or less than 90% (as indicated) 
of the density at the center or undisturbed region of the plasna. 
2. Plane and Rectangular Configurations 
ae Rectangular cavity 
Because it is the most complex geometry we can treat in 
rectangular coordinates, we begin with a discussion of a rectanguiar 
cavity of dimensions X95 Yo, and Zo. The limitations of the develop= 
ment of the diffusion tensor prevent the application of the followe 
ing solutions to magnetic fields, and thus the scalar coefficient D, 
will be used. Of the possible variations of the solution (such as 
equations (3-11) and (3=-13)), we can immediately rule out those con- 
taining the exponentials because of the requirement of zero density 
at the walls. The boundary conditions also limit the solution to 


the first positive segments of the sine functions (and exclude the 
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cosine functions) if the origin is taken at one corner of the cavity. 
If the origin were selected at the center of the cavity, the cosine 
functions would be used and the sine functions excluded. With the 
origin at one corner, and the sides along the positive coordinate axes, 


we obtain the following solution for the steady state (k, = 0) plasma; 


7(4,4,2) = Yo sin{Z x) alt 4) wn(Z2) (3-52) 
with the diffusion length given by 


pe ‘ 2 
ee Cl ca) ee 


(3-52a) 

The three components of the particle current density at any point 
are easily obtained, using equations (2-2) 

= -me Ds E coal Ex) an(By) ome( Ee) (383) 

with similar equations for the other components. From these equations 
it is seen that the particle current density is zero at the center of 
the cavity, is directed outward from the center to the walls, and is 
a maximum at the walls where, for the data of (3-50) and a 5 centi- 
meter wall-to-wall separation, the magnitude of the current density 
at the center of a wall would be approximately 2.7 X 1916 particles/cm@ 
-sec. It is of interest to note that the steady-state particle current 
densities at the walls decrease in inverse proportion to a linear in- 
crease in cavity dimensions. This decrease is a consequence of the 
decrease in the gradient of the current density. From equation (3=52a) 
one can see that the rate of ionization, VY , necessary to replace the 
particle losses to the walls is inversely proportional to the square 
of the dimensions (and, incidentally, is primarily determined by the 


smallest dimension of 
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the cavity as a consequence). The explanation for the inverse 
dependence of the rate of ionization per unit volume is easily seen 

by equating the number of particles destroyed at the walls to the 

total number created within the cavity. The latter is proportional 

to the volume and to Y . The number of particles carriec to the walls 
however, is proportional to the surface of the walls and to the current 
density. The latter is inversely proportional to the linear dimen= 
sion, as previously mentionec, and consequently the number of parti«- 
cles drawn off is proportional only to the first power of the Linear 
dimension of the cavity. Consequently, the inverse dependent of 

on the linear dimensims is expected. The fact that the number of 
particles destroyed at the walls is directly proportional to the 

linear dimensions of the cavity (and therefore, so are the eneray 
losses) is alsc worth noting. 


The space charge field, according to equation (2=36) is 





Es = (Rr \neablsr Lad tbeagye] os 


or, again using the data from (3-50), is approximately =270 volts/n 
(directed away from the center of the cavity) at the point (.25x,, 
o5Yos e5Z9)- The cavity dimensions are assumed to be Sen. Again, 
an inverse relaticnship to the wall-to-wall separation exists. The 


Space charge potential (equation (2-;1}) is 


Vc Dabs si bien  fn§ oinlE.»\ an By) ainl (2)} (3-51. a) 


or “eens =-1.53 volts at the point (25x), Yo: 252.) for 

a5 om. cavity. The potential is with respect to the center of the 

cavity. The petential predicted by equation (3-5) becomes singular 
Ke) 





at the walls, and this feature of the solution indicates a limitation 
of these equations: they are not applicable beyonce the sheath edge. 
Using the current equations (3=53) we can compute the energy 


lost at the walls. From equation (2-101) we obtain 


Eun = 2n,D,(eVir EL, , wn E,x\onlE y)dydy + 


Ee Volts) Th e\kde + B E(x 3 2 plan edge 


where the factor 2 takes into eccount the surfaces adjoining the corner 


(355) 


opposite the origin. This equation reduces to 


Sn. De(eVit KE 
n a u ) ete | Xa2e 4 Xoty | 


Wo te 


_ BaDl XB) [ te Sa 4 x,% 4 ae (3=55a) 


The heat eerste per unit area is ini by dividing each integral 





wall = 


n (3<55) by the appropriate area, obtaining for the three sides of 


the cavity: 
dn Ds (eV; + =) 4, Ds(eVi+ Kr dine Ds (eV, 4 EE) 


.* TI Xe > M4 ) TI Zo ga”? 
Again using 5 cm. dimensions and the data from (3-50), we see that 


Cua] j iS approximately een 1017 eV/sec 27.8 joules/sec. From 
equation (3=55a) it is seen again that the total heat lost per secnd 
is proportional to the linear dimensions, and that doubling the 
length of all the sides, for instance, will double the rate of heat 
loss. 7 
One very important fact about ns characteristics may be seen 
from the preceeding equations, and in particular from equation (3=52a). 
That is that the effect of a small dimension on the characteristics 


ig very much greater than that of large dirensions, so that a long, 


narrow plasma has, in the middle region of the longer dimension, 


7 








characteristics that are almost completely dependent on the shorter 
dimension and are relatively independent of the longer dimension. Thus, 
a plasma whose z-dimensicn is very much greater than its dimensions 
in other directions can often be treated quite properly as "infinite" 
in the zedirectien, with corresponding simplifications of the equa- 
tions describing the plasma. 

The zones cf influence of the rectangular surfaces are indicated 
in Figure 3-1. 

If ionization ceases, )= 0, and the plasma decays exponentially 
in accordance with equation (3-5), with ky given by (3=7a), (3-12), 
and (3=-52a) as 


_ 
ey > a : (%, ‘+ (3,) ee: 


¢ 


Thus the plasma will decay to 1/e of its initial value in time. 


| 
D.\e)* He. é +(2\5 a (3=57) 


which, for a cavity of 5 cm dimensions, is about 20 Msec. This 
decay time is inversely proportional to the square of the linear 
dimensions. 

If the iritial distribution is not the steady-state distrim-= 
tior, but instead some function f(x,y,z), then the coefficients of 
equation (3-11) must be obtained in the normal manner for a Fourier 
series, and a series of the form (3=10) must be developed. For example, 
if the initial distribution is a costant density n, throughout 
(except discontinuously zero at the surfaces), the particle distribution 


function will be 
22 Abe DE 7 
20t/)nx : (2mm +i) n: - | Lin 
Aly 4,2;t) = Ef on Ie PNG ligee eee er | 
118 (3-56) 








with 


ues 
(k,), =| ai [emcdn |" fe (3=58a) 


It is to be noted that the higher order terms cecay at a faster rate, 

so that the plasma approaches the steady-state distribution as it decays. 
Figure 3-2 shows the decay of such a plasma at three different tines. 
Other features of this solution are that the rate of decay is in- 
versely proportional to the square of the linear dimensions so that 

a plasma ina large cavity decays more slowly than one in a small one. 
Further, we again see ee the shortest dimensions control the decay 
rate, so that dimensions which are very much larger than the smallest 
dimension may again be treated as infinite. For instance, for a 


length which is ten times: that of the others, the error in ignoring 


it in equation (3-57) is less than one percent. 


b. Infinite tube, rectangular cross=section 
By allowing Z2, tc go to infinity and setting f= a 
in the equations cf the previous section, we obtain the equations for 


an infinitely long rectangular geometry. Thus we have for the steady 


state 

(4,4) = Noam(H, x)an(E 4) ie 
with 

| om 

+ kt ae, cy (3-592) 


The current density, potential and other equations are similarly ob- 


tained and will not be stated here. 
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c. Infinite Plane. 

_ By again allowing one dimension to become infinite, for 
example, yo, and by replacing y by y,/2, we obtain the equations for 
the simplest geometry, two parallel planes separated by a distance x); 
or a Single plane in contact with an infinite plasma in which 

x, = 5 Or (3-60) 
where we know the rate of ionization from other considerations, or 
where x,/2 is the distanceto the undisturbed region of the plasma. 
The restriction that there be no magnetic field can now be dropped if 
the only density gradient is along the coordinate lines perpendicular 
to the planes and if the field is parallel to these planes. If the 
field is perpendicular te the piane (so that it does not affect diffusion 
perpendicular to the plane), or if there is no field, then Dry in 
equation (3-60) and in the equations to follow can be replaced by its 
no=field limit, Ds. 


The particle density distribution for the steady state is 


N= Me oun(35X) , (3-61) 
with 
> r= 
i? = ae) 610 


The other equations are derived from those describing the rectangular 
cavity by replacing y and z by Yo/2 and Zaye respectively, by letting 


Yo and 2) become infinite, and by replacing Ds by Dm) 


3. Cylindrical Configurations 
ae Finite Cylinder, interior 


For a cylinder of radius R and length 2,, we can immediately 


Simplify Equation (3-16) for the steady state by considering the 
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functions and the boundary conditions. The Neumann functions are 
singular at the origin, and so their coefficients, C3 Oy? must be 
zero. Further, if the origin is taken at the center of the Ease of 
the cylinder, the coefficients of the cosine functions, Cans must be 
zero also, because the density is zero at z = O (using boundary 
condition (2-99), Then, as in the case of the rectangular configura- 
tions, the requirements that n always be positive (or zero) and the 


boundary condition at z = zo requires us to set all the sine 
coefficients to zero except Cb p> and Ky becomes z - Because we 
have discussed no forces which would tend to establish an assymetric 
angular distribution, and indeed in the case of the derivation of 
the diffusion tensor ir a magnetic field have specified that no 

such assymetry exist, we set m= OQ. Thus we are limited to the 
zero order Bessel function, and because of the boundary conditions 


must select the first positive segment of that function. Thus we ob- 


tain the well-known distribution function 


nls2) = Me Tol r)ouwl Ez), (3-62) 
where 

Re ret ae (BY + (), (0 
and 

ye (*R ) Dm, + (% g.) Drie (3-62b) 


The magnetic field, if any, is parallel to the cylinder axis. 
The 10n current density components are obtained by using 
Equation (2-2): 
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(3-53) 


Thus, at the center of the cylinder (r = 0, z2 = 4 2,), the current 
density is zero, and within the remainder of the cylinder i and ey 
are both directed from the center to the wine: © is a maximum 
at the ends of the evident but l recches a maximum (see Figure 


3-l) when 


incopy at r = 0.7SR- This maximum occurs at a radius less than R 
because the additional current supplied by ionization beyond r = 0.75R 
is not sufficient to counteract the effect of the increasing cylin- 
drical volume into which the current eeeremiae, Banc fen the rectangu=- 
lar cavity, an increase in cylinder size is accompanied by a decrease 
in ionization rate and.current densities in the steady state. 


The space charge: field, as obtained from (2-35) and (2-36) is: 


; . —p.. ae 
ry “Ee inet 2) ow 3 ae a é °) 
| : <r) (3~61:) 





eee, Do = De We sA(L2), 
The zones of influence of the cylinder sume are indicated 
in Figure 3-3. 
One of the characteristics of this plasma which is of interest 
Pomthemtaet that the ehate of the beenay state configuration of ion 
density is dependent only upon the cylinder dimensions if the strength 


of the magnetic field is less than that at which the "pinch" effect 
5h 





begins to appear. However, the ion current density in the radial 
direction (perpendicular to the field) is proportional to the inverse 
square of the magnetic field (at high B) through its dependence on 
Dm, 

If the steady state plasma begins to decay, the distribution 
is given by the equations 


N\é, 2.¥\2 
- Kz) Om =x } Dw 
Wr, 2, t) = ee 2.) 3(7r) 2 Kz) i sll R rte (3465) 


If, however, a decaying plasma with initial distribution F(r,|f ,z) 

is to be described, it is necessary to develop a series from Equation 
(3-16) in the same manner as used to develop a Fourier series, obtain= 
ing equations similar to (3-58). 

be Infinite Cylinder, interior 
In a manner similar to that of the rectangular problem 

it is possible to treat a very long cylinder as if it were effectively 
of infinite length. Employing identical procedures, we obtain for 


the particle density distribution function of a "long" cylinder of 


radius Rs 
oo TA) (3-66) 
where | | 
ri : (=e) ond Ys (2*) Dm, , (3-662) 


The particle current density is given by 


4 
[! = + Dm, nT, (Zr), aie! (3-67) 


r 


the field by 


[(-42B°){Dm, -D- 2.4 +l) 


eS ye -yuTt nO. (3~68) 
ae Me R TEs) ( | 








the space charge potential by 


os oo nae Dir 
V = ey, (7 — J (3-69) 


and the heat liberated at the walls per unit cylinder length by 
T- tet 
(3-70) 
1.85 7e Der, (eVia iy) 


It should be noted that the latter is independent of the cylinder 


{! 


Evol / tanith 


{> 


dimensions, anc using the data of (3-50), is approximately 1.32 x10? 
B7¢ eV/m-sec = 2.11 B"* joule/m=sec (with B in webers/m°) 

The decay problem is treated in the same manner as for the finite 
cylinder, and we obtain for aecey from the steady state 

mrt) = Ne 3 aot -(R = Dt ( 3=71) 

Figure 3-l, is a plot of the radial current density as a function 
of radial position for a steady plasma. Figure 3-5 is a graph of 
the space charge potential versus radial position for an infinite 
cylinder. The fact that the space charge field and potential become 
infinite at the wall is a consequence of the assumption that n,77 
= QO, As in the rectangular case, these equations are not valid within 


the sheath. 


Geb Ent inice Cylinder in Infinite Plasma 


Using our techniques, there is no exact Steady-state 
solution to the continuity equation (3-1) or (3-2) for a plasma exterior 
to a cylinder of finite length because we cannot combine the solutions 
(3-16) in a manner such as to meet the boundary conditions at the ends 
of the cylinder. However, for the infinite cylinder of radius R 
it is possible to solve for the distribution function. Again ve limit 


ourselves to those problems in which no assymetric forces are acting 
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(and thus to axially symmetric plasmas), but retain the Neumann 
functions because the origin is no longer within the limits of r. We 
obtain from (3-16), for the steady state, 
tate a 
u(r) = A fx.( y+ CNAE, vr) . \xeree (3-72) 
where A is a Bteereer toh constant which we shall determine later, 
R is the cylinder radius, R!’ is that radius at which Mr) .1s 3 maximum , 


and 


TJ be AR) (3-72a) 


; : NE i 


if the boundary conditicn (290+ (nNy911 = 0) applies, and is a more 
complicated function if equation (2-71) applies. The latter function 
will be developed in the section on probe theory. R}' the radius at 
which equation (3-72) is a maximum, is found by setting the first 


derivative to zeros i.e., it is the solution of 


(E.R) = on GER (36720) 


The normalization constant, A, can be determined if the density n(r,) 


Boa radius r, is known. Then 


aia (3-72c) 


T(J) + ONE Dey ch 


The radial Li enone is 


[= Do Aix, 4 5 IZ r+ cn (Vz. “f ou 


-_ 


d. Coaxial Cylinders 





The equations of the previous section also describe the 
particle anc current density relations which exist between the surfaces 
of coaxial cylinders of infinite length with inner radius R and outer 
radius R!' , Because of the boundary conditions, Y is no longer an 


independent variable, but is an eigenvalue in the steady state plasma. 
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Using the boundary condition of zero density at the walls, the value 


of fz is seen from equation (3-72) to be that value for which 


T(E, 8) + CNR) =O (3-7h) 


h. Spherical Geometries 





Distribution functions for decaying plasmas with arbitrary 
initial distributions and for steady-state plasmas in Spherical geo~ 
netries can be obtained by series development of equations (3-17). 
Equations for the spherically symmetric plasma in the steady state 
(with ionization occuring) are very similar to (3-66) and (3-72) with 
spherical Bessel functions (equation (3-17c)) replacing the cylindri- 
cal Bessel functions in equations describing spherical geometries corres= 
ponding to the cylindrical geometries treated. In addition, it 
is of interest to describe the steady state solution for a steady or 
very slowly decaying plasma in which there is no ionization occuring, 
and in which a sphere has been placed. We treat the latter problem 
rst. 

a. Sphere in Homogeneous Plasmasno ionization occuring 
If we have a spherically symmetric infinite plasma 


with an undisturbed particle density, n 


0? in mich Y= =n = © ; 


¢ 
with neither magnetic nor electric fields imposed, we have from 
equations (3-17) (with kg =m = O) the following density distribution 


around a sphere of radius K 


A(t)= No ~ — re (3=75) 
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We have employed the boundary conditims that n(r) be zero at the 
sphere and equal n, at large distances. The zones of influence 
for such a sphere are 

BO, = “ere 


90% - Y<IOR (3-76) 


) 


and the radial current density is 


ft. - De %eK G77) 





nD 
This current density is a maximum (of magnitude —2 =) at the surface 





of the sphere and approaches zero as 1/ré at large re The current is 
directed toward the sphere and represents a flow of particles from a 
source at infinity to a sink at the sphere. 

The heat liberated at the surface of the sphere is 


kha 
s 
as 


Ealune ~ E NR ev, De (eV; + a (3=78) 


It is proportional to the square of the radius of the sphere. 
Using the equations thus stated, it should be possible to calculate 
the particle density distribution which a small spherical surface in 
a large, slowly-decaying plasma would establish about itself. 
be Spherical Sheli-internal plasma 

From equations (3-17) it is also possible to obtain 
easily the steady state relations for a plasma bounded by a sphere of 
radius R centered within a spherical shell of radius R'. whe shall 
again assume Spherical symmetry, thus limiting ourselves to the 
Spherical Bessei functions of zero order. By employing the follow 


ing relations for these functions 


a 
cooler) (3-17) 
¥lo Ck) - 7 “er 4 
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and doing some manipulating, we obtain 





Ti ee ie | 


Fe Cie pie| pe i eeney, ( 3-80) 
wnere 
ra 
yee! ax) i 
Be _ /A* aod (—— G3 80a) 


and for which derivation the condition of vanishing density at the sur- 
faces was employed. ©, is an arbitrary normalization constant which 
can be determined if the density is known at some radius r, within 


the plasma. 


Then, 
| 2 f TieaR5 : 
A, = Gig LAPR Lo ( (3-80b) 
and 


eo = “er, CSc | ak =F a B| (3=80c ) 


/ 


The radius of maximum density, R't, can be found by setting the first 


derivative of (3-80) to zero, which yields the following equation for 


cE cox | TL =k i poof ES | (3-81) 


For a sphere within an infinite plasma in which ionization is occuring, 


Rit 


this is the radius at which the effect of the spherical surface be-= 
comes Zero, 


The current density within this plasma is 


= -71r Ds se [A je i 


at jan |- suns a(r-P 
x AB CR), rare te. po et ee) 





(3-82) 
It is to be noted that it has two components. The first is directed 


toward both surfaces from a radius within the plasma, and the second 
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is directed radially outward from the center, but is zero at both 


surfaces. 


For a spherical shell of radius R! with no internal sphere, one 


obtains 
Mok! uc 
Mr) = Te sa ; (3-83) 
with 
es ee! 
Ke Fe (3-83a) 


This equation was obtained from equation (3-80) by setting R =O, and 


by evaluation of C, by L'Hospitalts rule 


hum - CG, am ton ( 35) , Aste G z coa(2E ) =¢ 


f-O —— | =o (3-835) 
where n, is the density at the center. Figure 3-6 is a graph of this 
function. 


5. Prolate Spheroidal Geometries (the wire or probe) 
The influence of a prolate ellipsoid upon a plasma is now 

considered (see Appendix I for transformation equations and figure). 
In its limit, when the minor axis is small and the major axis approaches 
the length d, the effects of such a body can be used as ar approxi- 
mation of the effects of a finite wire or narrow probe of length d. 

Flammer (6| shows that the wave equation 

Wes k*) N =O, (3-8) 


equation (3-7), expands in prolate spheroidal coordinates tos 


5 (iy) 5y + 53 SNE + eat see tA TIN = 0, (3-85) 
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in which we have set 
Oe 2 kd, (3=85a) 
and that, by the methed of separation of variables, solutions to 


(3-85) may be obtained in the form of the Lamé products 


Nena = Somn G7) Reon (6 Sen (Og (3-86) 


The functions Snn(es7] ) and Rmn (c,? ) satisfy the ordinary differential 


jn ld San (67 + Amancty® - a Sal) =O (3.86a) 


FLT Ne Ral Gl Laen CF + 7 | Rin (6 $)=O (3-860) 


in which the separation constants ae and fy are the same in both 
equations. It is a combination of such functions (of the first and 
second kinds) that is given in Equation (3-18) as the general solu- 


tion to (3=8),) 
a. Prolate Eilipsoid in stationary infinite plasma 


Before proceeding further with the prolate wave func= 
tions, one special case, analogous to that of the sphere in the 
infinite stationary plasma, should be treated. By imposing the 
boundary condition that at large 5 the density, n, of a stationary, 
non=1onizing plasma be independent of \? and 4 (so that SH = 0 


2 
and ar = 0), we can reduce equation (3-85) to 


te lsat *o (3-87) 


which integrates tos 


: i 
él 





where no is the uniform density at infinity, and the density at the 


surface . is zero. The current density is then: 


I ™ cot “'(6,)(3" om 1) 


To investigate the properties of this solution, we shall use 


(3-89) 


the identity 


cS) = FE dnlhtl) - GLal§-I) (3-90) 


by which Equation (3-88) becomes: 


; — ghnl$St!) — £.OnGg -1) : 
1 =e EOn(Sot1) -— dnl S, =!) ee 


at very large distances (i.e, § >> 1), In (F +1) = 1n(F 1) “a 
and equation (3-91) can be expanded to 
Se 

Ax Heol | a ) (3-92) 
which is identical to equation (3-75) for spherical geometry, a 
result which reflects the fact that the surface "looks like" a sphere 
at large distances. In the vicinity of the surface, i.e., for omer 
we let 

$ = o +6 (3-93) 


(with €such that the mean free path, d » is much less than € ), 


and expand, obtaining 


3  £n§.r€ +!) - saLte )| 
" 7 kn(Forl) —nlt, -1) ae 


é5 





-n{- Bnlgct) + dnll taj) -Inl§e-l) - Sal ¥& ) 
_ _ (3-9hb) 


ee +1) Ae I) 


E€ 
M1 +55) - L(+ 35 &9ic) 


=Yle 
bnGtt) -4uGeo-!) 


s 


By neglecting second order terms in € and using the expansion 


x2 
Gaxe!\ = nea Sy bore, (3-95) 
equation (3=9lic) can be reduced to 
ee 2 
at oc J 
=} |, ewe 3-96) 
NA No hn ie 7 No e bn er 
an Oo” : 


Thus, for the density at the surface, € = 0 and therefore n=0. Equa- 
tion (3-96) shows that the density is a linear function of the 
distance from the surface for short distances. 

At large distances, the current density can be seen from Equa~ 
tion (3-92) to be an inverse square function of distance (as it is 


for a sphere): 


re Ue. 
is = - DV = afk (te) at —_ 


At short distances, we obtain from (3-96) 








= : 
v7-_— (3-98) 
: “eye Gea ES 

which is a minimum at the center of the SS ELL) f 0) and approaches 


infinity at the ends C « £1). This occurs because the area at the 


ends is zero, and thus any current flowing towards them constitutes 


eo 





an ‘infinite! current density. 


The current density at the surface is found from (3-98) to be 
~Zne Vs 
wa a a 
f° (G2 -)CB-79*) BS (3-99) 


which again becomes infinite at the ends and is a minimum at the 
center. We can obtain an estimate of the heating at the surface 
per unit length if we approximate the ellipsoid with a cylinder of 
radius r = (f - 1F and approximate the current density with that at 


the center (7 = QO). 
L-T- 
Baar 4M n(e Ve + a) 
Uumat Lernrath Sa 
d Zo On 
To find the behavior at small radii, ve let ¢ =1+é€, obtaining: 
Unne (eV; + =) 
E Set iste, Lh. a ee 


wy ~tiVMe (eV, +i 


(3-100) 





(3-101) 


Eo 
which approaches zero as S. approaches Zero. 
b. Prolate Spheroid_in steady state plasma = ionization 
eccuring 


Given an infinite plasma of homogeneous density n, at 


O 
large distances and a prolate spheroid whose surface is the surface 
yy one can simplify equation (3-18) considerably by consideration 
of boundary conditions. First, the requirement that the wave function 
be finite at Y | = + 1] confines the solution to angle functions of 

the first kind since the angle functions of the second kind are singular 
at these points. If we look for those solutions for which the der- 


sity, n,is uniform at large distances for all , then N is a function of 


4 and ° alone, andm=0O, Further, the prolate angle functions 
E7 





reduce to the associated Legendre functions when c vanishes; i.e., 
as the sphercidal insulator becomes spherical in shape. Therefore, 


a 


we can determine the subscript “n’ in Equation (3-16) by studying the 
Legendre function, Pal 7] ) because our solution must satisfy the Limit- 
ing case of a hemispherical body. Only Pa 7 ) is positive for all | 
in the range -1l> i }. Because n must be positive for all a the 


subscript must be 0: Consequently, we obtain for the density distri- 


bution 
M = 1, Soe G )|Ree (c,2) ‘ BR. (<, a], ( 3-102) 


where 


oe Pies (3~102a) 


Tf we employ bouncary condition (2-99) (N,.3; = 9); 


= anon) 

b= - Rows) (3~102b) 
It must be remarked that this solution does not satisfy the require-~ 
ment that the density be independent of 1 at large distances. However, 
for small values of c, and therefore for either very small or for nearly 
spherical bodies the 1) dependence is negligible. The deviation is 
Only 1% for c = 0.5, but increases to practically 95% for c = 520. 

c. Prolate spheroid and infinite wall 
Let us consider a semi-infinite plasma of uniform density 

no at large distances, bounded only by a single infinite plane which, in 
prolate spheroidal coordinates, is 4g O. Purther, consider an 
insulator (probe) whose surface is defined by S=%, projecting into the 


plasma as indicated in Figure (3-7). The solution to the continuity 
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Figure 3-7  Prolate Spherdid ond Infinite Udall 





equation is developed in the sane manner as in the previous section, 
except that in this instance the appropriate Legendre function is P97 ) 
because it has the property of being positive for all values of 

in the range 0 to +1 and further of being Zero at 4 = O. Thus the 


density is 





(1) cana (2) 
Wie 715 as Cn) LR (Go) a 3 Ko (<, 5] (3-103) 
where, again, 
kd dAlv a 
and if N11 = 0; 
(1) 
B=- Ror ge (3+103b) 


K a (<, $0) 


The field and potential equations along the axis of the probe 


[ie., for ns 1) ares 


ee [Rov (o$) + BRO, $) 


a (3-104) 
=| da te = i a 
eA Mer FRO Ce) BROCE] 
where the prime indicates the first derivative, and 
. ue Dee. Ny 
Mabie a ad Haag. a (3-105) 


= 2 Ded. La} Railef.) + B Rei Gi)? 


Mp CRAGS) + BRIE?) 
6. Oblate Sphercidal Geometries (circular lamina) 

In general, the solutions to the continuity equation in 
oblate spheroidal coordinates are the same as for the prolate with 
-ic replacing c and +10 replacing ° in the equations. Consequently 
these will not be further discussed, except for the problem of the 


circular lamina in an infinite, stationary, non-ionizing plasma. 
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a. Lamina in stationary plasma 
To study the influence in an infinite plasma of a 


circular lamina, the effect of an oblate ellipsoid is first discussed 
amd then the limit as the semi-minor axis approaches zero. 

Proceeding as for the prolate case, by imposing the condition that 
at large ° the plasma density be independent of ¢f and y » and by sub- 


stituting +i for % we obtain from equation (3-85) 


~ 
fers )3E | =e (3-106) 


Integrating twice, we obtain 


y 2G ton (2) +C, 
Letting . N =NMs 5 we get 
y= Cp tan CO) Co = ot eee (3-108) 
Using aS the boundary condition 
M5.) = Ce tax (3) tC (3-109) 
subtracting this from equation (3-108) yields 
Yo 
- LSTi- toa &) o (3=110) 
=| 
Alo Tom ($2) (3-111) 


ie Cz — 


1.$ 711-0 (3) 
so that, for an oblate ellipsoid in an infinite plasma, 
taw'(#) = ton'C§) 
A= 20} : (3-112) 


S71 = dam'( So) 
a 





When the semieminor axis is zero, 5 - 0, and the figure is a 
o 
circular lamina of radius ry, = c. Equation (3-125) becomes 


No tare (3) 


The zone of influence of this mina, within which n =3Nno, is given 


by 
tam, (2) 2 nen = 0.785, (3~11h) 


eC... f= [ (3-11ha) 
which is an oblate ellipsoid with semismajor axes r = joo and 
semieminor axiS Z = To. 


The current density in this plasma ( with c ’Y; ) will be: 


* -Dynot JER 5a 2 Pete (3-115) 
"\ LSnG (s2meVge+l) 
This current density is a maximun at the surface of the lamina, 


decreases as to zero towards infinity, and because the area of 


\ 
cs 
the "side" of the lamina is zero, is infinite at the sides (7 =0) of 
the "surface", If we assume the current density at the center of the 
surface to be an approximation to the surface current density across 


the entire surface, then the recombination energy liberated per unit 


surface will be approximately 


TT Da Not leT- 
E apie XY ee (eV: F Lay (3=116) 


which, as for the sphere, is proportional to the radius. 
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IV 


APPROXIMATE SOLUTIONS FOR INSULATED PROBES 


When an insulated probe is inserted into a dense plasma--that 
is, a plasma in which the mean free path is equal to or less than 
the probe dimensions~--Langmuir's probe theory 8) may not be appli- 
cable. Particle losses at the probe surfaces may cause appreciable 
effects on the particle density distribution, and hence on the probe 
potential. Under these conditions, « of equation (2-71) can be 
appreciably large and consequently so may the effect of the in= 
Sulator on the steady-state distribution. The approximation describ- 
ed in Section III=A=2 has been employed to calculate this effect 
for a floating probe in a partially ionized gas. Two different 
geometries are considered. The first is a prolate llipsoid 
inserted radially into a cylindrical discharge, with the sensing 
element on its tip. The other is a cylindrical probe aligned parallel 
to the discharge axis. For the latter, we are able to account for 


the effects of a longitudinal magnetic field. 


i 





A. The Prolate Probe 


We do not have an exact solution of the continuity equation 
(2~55) for a cylindrical or ellipsoidal probe projecting from a 
cylinder wall into a plasma. However, we have seen that there 
is a solution for the density distribution in a long, narrow 


plasma (equation (3-66), In the absence of a magnetic field, it is 


N=No Te(&r), (et) 
in which 
. — 2405 (h-1a) 
aan KR. 


| 
Fiqure bbe] Holate Drabe I Cubndrical Discharge 
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Further, we have obtained a solution for a prolate spheroid attached 


to an infinite wall (equation (3=102)) 


N= Mo! Soil, ”) [Roe £)+ BRar Ge, (11-2) 
Based on the assumptions outlined in Section III-A=-2, we shall 
approximate the effect of a narrow cylindrical or prolate ellipsoidal 
probe by joiring the two models described above as indicated in 
Figure h-1. The infinite wall is tangent to the cylinder, and the 
z-axis of the probe (7 = 1) lies along the diameter of the cylinder. 
We assume that the presence of the probe does not significantly alter 


the ionization or diffusion rates, so that equation (3-103a) becomes 


dfvy ~ 2.42 
Cz LE as te (hi 2a) 
where wf, the length of the probe (if the minor axis is much 


¢ 
shorter than the major axis) and R is the radius of the cylinder. 


No effort will be made to describe the density throughout the 
cylinder. We shall merely attempt to predict the density distri- 
bution along the 7 = 1 axis of the probe (i.e., along the diameter 
of the cylinder). It is this density distribution which is res- 
ponsible for the potential distribution sampled by a sensing element 
at the tip of the probe. In accordance with the procedure outlined 
in Section III-A=2, we readily obtain the following condition for 
determining the position along this axis of the dividing surface be= 


tween the two models 


Stlacane “Hh, 
“ 'Probe Cylinder (hi=-3) 


Equation (l-3) was solved graphically because sufficient published 
data on the derivatives of the spheroidal wave functions was not 


available. 
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_ Because of this paucity of information on the spheroidal wave 
functions, it was also necessary to limit the study of this problem 
to the limiting case where d is so small that the boundary con= 


dition (2-99) (n = QO) may be employed. For this problem, B of 


wall 
equation (li-2) is readily obtained: 


(1) 
6 = = 
Ro (Ga 


The normalization constant, Nos of equation (l-2) was determined 


(h-l) 


in relation to n, of equation (l=1) by the graphical solution of 
equation (l-3). 

A study of the solutions, and of the physical properties of our 
model, reveals that this model can be justified as an approximation 
to the particle distribution in the vicinity of the probe for probe 
penetrations varying from zero to about half the cylinder diameter. 
For short probe lengths, particle losses to the adjacent regions of 
the infinite wall (which, for the short probe, is a reasonable approxi- 
mation of the adjacent cylinder wall) is the predominant factor in 
establishing the particle distribution in the vicinity of the probe. 
For instance, for a probe length of 0.2R, the parameter c of equations 
(4-2) is 0.48. From Figure II-1 (Appendix II) we see that the pro- 
late angle function sQ) is very nearly that for c = 0, which 
corresponds to zero probe penetration. Consequently, the effect of 
the short probe is very slight in comparison to the effect of the 
wall, as we expect from consideration of the surface which would be 
the cause of most of the particle losses from the region around the 
probe. However, for deeper penetration, the angle function changes 


appreciably, indicating that the surface of the probe exerts a more 
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dominant influence. Again, this is expected since particles in the 
vicinity of the tip of the probe are more likely to diffuse to some 
portion of the probe than to succeed in reaching the wall. Thus, 
even though the infinite wall is no longer a good approximation of 
the cylinder wall for the deeper penetration, the increased probe 
surface losses justify using the resulting particle density distribu- 
tion functions to describe the plasma near the end of the probe. 
When the probe penetrates more than half-way through the plasma, 
the influence of the opposite walls of the cylinder on the particle 
distribution near the end of the probe become increasingly more important, 
and our model breaks down. 

Conversion from prolate coordinates to probe dimensions is accomplished 


with the following equations. 
A 
probe length: ff = 3 S zt (x5) 
probe radius ( = seni=minor axis): ne Aes Tp ek fe ae 21-6) 
P 20 oI Py ; 
Re he is the equation of the probe surface, 


Flammer |6| was the source of the information on the spheroidal 
wave functions. Data not calculated or graphed in this reference is 
presented in Appendix II. 

The graphs on the following pages were computed using the proce- 
dures otitlined in this section and subsidiary equations developed in 
preceeding sections. 

Figure l=2 is a graph of density versus distance for probes of 
Several radii but equal depths of penetration. Note that the magnitude 


of the dependence of the probe influence on the probe radius is 


(7 





relatively small, an order of magnitude change in probe radius 
causing an increase in the region of influence of only 50%. 

Figure h=-3 is a graph of the density distribution for various 
depths of penetration of a probe of radius Tp = O.02R (which, for a 
2.5 cm discharge, corresponds to a 0.5 mm probe). 

Figure hel; shows "effective probe length" vs. true probe length 
for the 0.02R probe. The effective length is based on the inter- 
section of the two models, and is approximately the point of the ori- 
ginal particle distribution that the probe samples. The curve marked 
rae ar, is the effective probe length when the probe does not 
disturb the plasma, and corresponds to the true probe length. The 
broken segment of the curve is an estimate for that region in which 
this model breaks down, as outlined above. 

Figure li-5 is a graph of the potential disturbance associated 
with various depths of probe penetration. This potential disturbance 
is defined as the potential drop between the point the plasma samples 
(the effective probe length) and the sheath edge. A sheath depth of 
O.O1R was assumed for this calculation. These potentials can be 
used with the effective length data to analyze experimental results 
or to predict them from the theory. 

Figure l-6 represents such a potential prediction. It is the 
prediction of the potential measured by such a probe assuming that it 
actually measures a difference between the potential at the maximun 
density 1egion of the plasma (Zero potential) and the potential at 
the sheath edge at the tip of the probe. 

A very significant prediction of the above theory is that the tse 
of an insulated probe in a plasma of short mean free path will cause 
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a definite asymmetry even in an originally symmetric plasma. Although 
these results are not applicable directly to a plasma in the presence 


of a magnetic field, some of this effect is still expected. 
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B. Cylindrical Probe Parallel to Axis 


Figure h=-7 is a top view of a possible insulated probe con- 
figuration for measurements in the presence of a magnetic field. The 
cylindrical probe (enlarged in the figure) is supported at its ends 
and long enough that the effect of the end supports are negligible 
in the region at the center of the probe. The sensing element (or 
elements) are located at the center of the probe on the common 
diameter of the probe and plasma. 

There is only one alignment of the probe for which an exact 
solutionof the continuity equation (2-57) is possible: when probe 
and cylinder are co-axial. As previously developed (equation 


(3-72)), the solution is 


n=7'| (ED - CN e) | SS aaa 


with lp the probe radius, R the cylinder radius, and n! the 


normalization constant. If A is not very much less than the probe 
dimensions, boundary condition (2-71) should be used at the probe 


surface. However, equation (2-99) (n = 0) is often still appli- 


wall 
cable at the cylinder wall, and was used for our calculations. 
Using equation (2-71), which is 
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Next, employing equation (2=99) at the outer cylinder, we obtain 


T.UE.R) = CNBR (etd) 


Equations (4=10) and (+11) are sufficient to determine both C and 


OR and they are most easily solved by successive approximation. 
ee ioe) 


Table h=1 lists the results for these parameters for various ™ between 
zero and one, using the following probe dimensions and mean free path: 


r 2=0,02R 
p 


r 


(again this corresponds to a 0.5 mm probe in aile> Cuacy linder), 


0.02R ] ' (li=-12) 


Cf 


Even when there is no recombination at the probe (c&% = 0), the rate 
of ionization for the steady state must be greater than when the probe 
is not present because the probe still occupies a finite w lume in 

the cylinder. 

When the probe is off=axis, it 1s necessary to employ the 
approximation described in Section III=A=-2. Equation (l=3) is still 
applicable for determining the separation of the zones along the 
common diameter (the only location in which we will be interested). 
In order to evaluate Pon, im the vicinity of the probe for varicus 
probe positions, we interpolate between the calculated value when the 
probe is coaxial and the value it would have when the probe is at 
the cylinder wall (and has negligible effect on the eigenvalue because 
of the low density there). This latter value we assume to be the 
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familiar 2.05/R. Assuming that this parameter would depend on the 
density of the plasma in the position in which it is placed, we 


obtain the following relationship 


oT zt) 
jae = “=| 2.405 4+ AY nn) (li=13) 


where 


Dy ~ Re, Voces —£.40S (l=13a) 


From the values thus calculated, a mean was calculated for each ™ 
(see Table h=-1), and this value was used in the computations of the 
density distribution in the vicinity of the probe. The use of a 

more precise value is not justified by the nature of the approximation. 
For the density distribition in the regions controlled by the cylinder 
walls, the value of a for the cylinder alone was used (2.),05/R) 
because the exact effect of the probe on this distribution could not 
be accurately estimated. The degree of error resulting from these 
approximations can be inferred from a comparison of the exact solution 
for the coaxial probe and the approximate solution shown in Figure 
ell. As expected from the approximation model, the maximum error is 
in the region where the two models are joined (where the positional 
error is 0.1R and the magnitude of the approximation is 12% greater 
than that of the exact solution). The magnitude of the error de-= 
creases towards either surface, but the percent error remains approxi~ 
mately constant. Cansequently, the potential determinations (which 
are based on a logarithm of a ratio of densities) will be reasonable 
approximations. 
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Several distribution curves have been obtained by using the 
procedure outlined above: (1) Figure h-8 is a graph of the exact 
solution for a 0.02R! co-axial cylinder for various < ; (2)Figure 
i=-9 is a graph of the exact sclution for co-axial cylinders of various 
radii at a given o{3 (3) Figure h-10 is the density distribution for 
various x for a 0.02R probe in an off-axis positions: and (i) Figure 
hell is for a given & and various probe positions for a 0.0! probe. 

Effective Probe Position vs. actual probe position is indicated 
in Figure h-12 and potential correction to be applied to probe readings 
in Figure l-133; we have assumed the sensing point to be on the side 
facing the center cf the plasma when the actual position is less 
than 4R from the wall. Potential curves are given in Figures heth 
and els. The former one is for one contact, oriented as indicated 
above. The latter is for two contacts, one on either side of the 
probe along the common diameter. Figure h-13, the potential correction 
curve,is just the potential difference between the sheath edge and the 
point of intersection of the zones of influence for various ¢ and 
probe positions. 

Figure h=-16 is a curve of heat Liberated at the probe for 
various values of probe radius when the probe is aligned with the 
center of the cylinder, and 217 is for various values of cc. It is 
assumed for these calculations that the fraction of particles which 
diffuse to the sheath reach the surface ami recombine, liberating 
(eV. + kTy ) of energy. 

The reSuits of the above theory indicate that the cylindrical 
probe should have considerably less effect on the particle distribu- 
tion functions than the prolate probe in a given cross~section of 


the plasma. This would be expected from the decreased surface 
88 





involved for the cylincrical geometry. Consequently, there is less 
assymetry produced by the probe, and the potential curves, particulerly 
with contacts on both sides of the cylindrical probe, should be more 
symmetric. However, the effect of measuring a potential at an 
"effective distance" causes some flattening of the predictec curve. 
Figure li-9 indicates that the density distribution in the vicinity 

of the probe is relatively insensitive to the probe dimensions. 
Consequently, the potential drop associated with this probe configura- 


tion is relatively constant for various probe radii. 
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TABLE l=1 CYLINDRICAL PROBE PARAMETERS=--EXACT VALUES OF 


v} andC ys, co for probe at center, and mean 
Dray 


values VSe OZ, 








Mean Values 





Probe at Center 





on | Qo 
R 

0 2eki 007 2h 007 
R R 





90 











tu i Peeen i 








ee Llepe gel eget TS LP lla ater lage alee ae Tie) 
7. 7.888088 
Pr] or 


Pit ert 
rw TA TT 


ime tT 
wt ty ee 


See 9 | RE pe ee eee etl alee eee ele eae ee ee eer emer aie eer ae ee ee 
fone | melee) SSeS e| ee oe mee elm ie lwo lol elaine laine oeingne hal tagey tele) el ey ee) ofa lea eel Tt lay }— I. 
ee oes) ama ape ele ope pepe: eee elo, ae eal ae eee ee ae eee eee LTT TAT Tere a 





Reet fale SIS iage tac) wee le seme ele ieee iat | en et etetet ffeil er ele) eases Seeley gee ee eae eee ern Cope Te ee tt Pe ee Fk yeni lee gt TR Pf 


S* fay s COC oe ee as gf) 
BECERRA ao pee ee a 
TTI ee itr Lit fifi Pi iitiiis ttt itt th) ho ee eee eae eae ee see ee aia) Lal ade eee ar eee ae fe) eee a Soleo ele alaie eae alelelarel_ (elelala,o) eel alanine egee ela) peg a lelel | fests eee are ECE ed 


Pea A P El totter lol gamete) |) alent repeater Ie 

cee eee is ieee ee Cea ee eee aie [ef pelea after tT Ee be a re eee GT OT Pt eet Takeo Ta erty rae A et folate Pe lee Sele elas) 
5 al mae eee ee ee aay eee ele yop fe eee ye | eek efete Paley tafe pe yy lol ha ee eee eee Te eee eerie eri) aires! le ee eee eee er ere a ele eee pe ae ef Rae a, fg fe egn gage oe Da 
een ea | fc eee ee | ge | eee ele alee le ieleyel lol fal | ae eee ol flee elie oagey TT ele lamel (oy | later ole] Se Belo alammpee | Coke t gg) Ter pe eb | Co ae eee Tee er a ee ae (| | 
P+ co RH HH re uane nace be RN est +] tt yo a et ee ve a aBE he co tt rtf Pip Att Ha Hat soaps ee Bera HH eee BBE A Set 

. AS Sat) ed j ei a 

Sas el eas) ele ee Ra Pee et Be Deelah bempeerh een en Mi ach ak Re ened Ce CT eae) aie ee ie ea) fer eek fea re eee Stet ty eae at a en | ope ae pale [ale [lenient lay 
& eee opaale ee melee eal al (ajepeiet selepepela | eeeiege | hele) pel alereare get at eee eee Ce Cee ee eee ft | red] orp 
L Se ee Pemp Bes  i e fe  f g  fagm p ge gn eg | 
Ce Pere rr eee [fe fa ya | a9 On eyo | eat (agg ef [oe PF ee ai] Pe 
poly i? le) (foie (ejamel_ ial Ge oe ae ee atte en TE yt ar Ono 8 aon ew fem QU | fmt 90 9 nie] of no (| cfm oa cnc cfg gen] P| 
Aen tip ieeel ieee leer | lel oatots(epe| | [lain otal, | tole CR tt a ee er Cee 
: : = a en Fe ft sf gE ge | 
ann gg ef ao |e 


© a meet ere (ne) | | nen eT Te ee riotapetel | 1 [oned 

Sa eb iereete | iain |=). | lo aieiee | S[e gla .e) 96m (oilers, lo( [oto . oe .ep6 6) 0,e eine noe ee ae oe eee ee eee ee 

Se ee 

BME Sie ee Cee CECE CECE HEREC EEE EERE ECE EEEEEE EERE EEE EEEEEEEE REECE EERE EERE ECE eee 

a ot A og p 
TTT Sey oe be TAP Ne eet ite) | eae gore | See) on) | ( (efor. areal] | ee eT eee ET eee 8 ee 
See nena ee eo eae one pedamimzmrchommapunjunsimmnanasharbentempon atpupndnisies iam bencanimivasiusiatican ethan aan lan ben/astsentp aaltepet anand ean ae a he al 
EERE EERE EERE EEE EERE EEE ERE ESET EEE eet ae SSS | oy 
Se 
Pea eae eel ae ope lle SoH Ee SSS 


Se es els) oe eee ee eee sol Ply fet feta, dey 4 
Pee en Saregama eee epee Ca es Tt oY aie ka ys 


‘eighteen ee So oot 
TTiIiilitAegiiliti tut 


aR De es ene aa aa eee etal | lope) yet peepee ys | ol (oleae he] 
es ee Sle ae ee eee eee ele sie elas) ayaa fale alela) (lela aie | | olelotaleyel lana] 
oh a a pa ee le eer le ep fale el) ole apt ee ala|agaa at ale |e pag neff] pm mee] 


BSE eee epee eee Pe 


eS EL) eae ee al ala oleate eee feel t alata lala kee ral lala palate lent let ae lel taal total") ot lafetl | ta] | 1). 1 tt) tt 4 
‘ + PtP | a a 4 eres IP po Sanne Sana a Dl cel a La) : Ct DN a cel eh ND Dace linea eee Rie eae cel ee ee eel eee liat eRe lee hace hee lreeles baal ve elahest 
; : ; ant 1 = Tr Sa 7 : TT 1 = 1 — . 4 7 — ’ " 7 P " Ce P ( “=? a " . , q 
= ea eae aml a ~«! _ — —— ——————— pee em — = —_— ll ee Le - = S a a ha - *_ a ane = i. a = s a ——_— = a 








Pet ae eee ie 6 eke ll ie) Lie io oer ke en tet or ae ea eee oe, i eel el EE Ve Bae et Te i eT iy io 
Ps TE es oa rt 
tie SIF ae 4 oS ri HE ti he Seo 

wii eee |) le ied PTT 

Pritt tie tee met) ere 


peti me er a ee re = = etd ae | ee det eae le 7EG SEeenEr, Peewee? 

lf >! Ct tl ae Y | oo | “Titi Lite TTT III rt dd ee I er a ee ee A a A ee 

4 — Lit tt + p+ + +H + +4 f+} Lee t+ 4+ a ttre tH We) espe Tastes pl TOY gia p asda) al ea toto ta Stott ye Pie ey eae faa Tata) fee ea TS Jin] Fea Spe te de Ale hd mete] Cee ee) | TT ae 
a 


) 
a 
Ld 
B 


ego peiome tes Ty 
ae ay ey hee) ae | Oe eee | ey ba eh glee eee) eee | epee eee) eee epee emcee | ae 


C 
s 
L 


‘ Pe ee ea 
Co oa ee ee ee a ee a ae ee, et 
i eee te ee ee ere ee ee ae ae a ee eer ae a 
et Ee le) eae mania, belek ol) alee) [elepel sale oe ales) Ty ee ae eee ee ofeielels) oe | ls ae oem ole, eee ae ee ee Pe ee ee I Tee Ee Te See eee eo eee ee 
eee eae ee ee es ee a ea dee) ea ale ale a ee eae ete late ee se siee epee mel wee eee |) oleate | a eee SSeS ee ee ee 
a 
a 


ae 
phate) | 
oe 
eet} _| 
Coe 





Hr 
fay. | 
ll re 
ar 





Peas 
bet 
tial 


a elaee ga peel ae ee alee eee see Cees Cree eile elelalceeieieielejelael lear Aeialaieye ls) 
ele ee eee ee lsiele ea sees mire srer Cet Lise aie) ale slelsleeielaie| ft lafalel [a | aereislergsts|: ie 
eta te ia el fae ete els ee Se ae Ee eee) ae ale ate eles See eee ae alate) ag 


RAE eee tee eee WY 


BBi GGGRE B83 GH ih 
pt ba eT eal ft ee) 
Pie eT 
| ae Pot ayy El 
Po 


eee ele sion Ree) fel) a eel tae TP ey a ean a eg 
Sed ee ee leer ee eee en ee eee eo 
eee eg em ee ee ek Bae hal ee oe eee ae ery ee eee eg lal 
hese efi, Solna eee Ce) ae ee ee ae Oe en aes 


$e ee n° eee i i fe] [a am a ge fe eff oem |e ke apg |p) 
Ses See eee eee seer aes Sieiealseiaies lee lleee | eieeisl wae leiseei ae 


hy 


5 
= 
H 
i 
7 
a 


wi 
CT es 
er 
= 


me ee see a EERE ee ee Ceaser eset fe ae ae 
ame ee tte ee eT ale ee eee eee ale Ear ey) eee flea | | fae faed PTT a 
Ce tee ee Se ee ee) ee ea 
See ee ISIS SEE eee eislels( ee sia ales Sis aan vee er aa eee a 
eee EE ees Se Ieee ae ieee ee) ea Se ey 
a ee ee ee ae eee ie 


OSS) feakates) oe ope aie Rls pie iene oie) [eel keg le eet sean ey ae 
ee ie ae eee alaleara ia ey ee eee TT | 
ee eee) ee) aes aa 


nae 


p td la aie ee alelae ee ley 
alert Bale eee tad A 
pet J Ce de oe TT 
it Eee 
Te ae leap) aad aie | ose eae aes [sll] ale fee] sh] 
Lier eT er To 
Pep SS uivgs i Te ‘aus flee] ag iia HH seam 
Pot i 
at 3 Ph meh] ope) [ { 
ae Ct TT ee f 
iT Re T  es H 
eee en een 


hl openers | | sey eee eee a 


F if SB Sl 3G URC SSRSR SSSS50Raea0 
A Lf camo a ey a as Fs Ge fe || pe a 
eee) see lelateeaiel eetaiel eas) 
pee feo) ls ete eh ea see [apf sen ee 
epee leays) 
4 Swe SEE GEE EEE Ere SCE Pee re re aamae 


rH 
aT 


— 


Pitre eee A ee 

PN St TP OE 
OG ye] SO a FR 2 Bean pepe NPD | FONT ||| eS et LT lala ere eee tT TT oe eee Pe Tp et tt 
Preto : ST Tritt tt) ee etal a oe ee E 

oe te ) ee ee el ee 
COE ee 
BEE ee 

’ SO a J 
Sitti ert Ne ST RIT TTT TET TT Pt CLL LCL 2 At ere Tol tae fet ates) tial tite] tapatre et 

BREESE EEA Roe caae PCC Ce ee CCL t eee ee eerie Eee 


i ee eae ee ee a ate eS ON 2.90) GSS 68588 [eel lela] leat [lege] Ly Se 


ee tape | [a areal a ere tae ala 
(hes) Sea ae eee eae aielal oie 


Lr EEE EEE EEN DD 
A te a AQ a tty oe an oe oo 4 ae a rt saaee gaseear-aaenes et De Ni I ae wi 
Pr Po pe ppp et ee fh py Ay id ‘oa 4680. ale eee Oe ae afale) pt TNA ON ON j 
CCE REESE ERE EE is ee a NE eR 
eee el eb dete ler eta bal OS ea a ee Ty fin ae Lot Se ll ee ete Ht oS 


Praeger ee ee 

















a 
~ : 
i } co B ed SB aEeears 
He : i Ro ee SBS S9 Seeee cere ooo suaEZeEES 
ch oe i — pont eecieedeenteehen a adm he) ao Ie Se Pk oe ee) i ooo 3 i. : 2B (Feat tt coed peek Se akc pomst age eo} 
. i - = 

= 


——— « 7 alia Sian a = a = 











a TT TPP PT ELS ITT ISI TT TT ee ewlenre ee lent eee] | fella | fee || 
Te et LP am tey ee fp une mpeg Poem en eel fern en] [ ery ahem Jeane] Pun fem yy Sate en eg eegeay | Peg ene ed Ley eye) een y ee |e ieee pene en eee lereley | lamer a en pal ee een ptf ttt Bee SEE eee Cer eaele aa eleie 

Yk] Lofmmemp de Reming enN emt |) mule) elope mT emNe | PemNEm No Neapem] P| Smee pss) — Jomtee) fem) enemy | gan ca bY] fen | Ofelia lease siete epatafel tr 
ele] Era aie Saree ar ease ian eee] | eae ee a eee el eel easieiel [opel ff laleielalafaiee fee Clair 





AD! Tt em TS TT eer a re eee es Po ! ; 
— — o er Y Q ? mac Triltlltt! Mililli ii eriri(y tea gy eT Pe we 
4 = otal Soon se Ey SS ++ Fae : pt } Stas Barents | ba | Lee Pe Hae seen fea a ezess PH A Fics [Vimy aml a YY a Ph Dall — (ol Tae ete) | ee ane | lee le! Loft) Sasa er ella vee we ebete? lot TS Slaaee 
ame eee ee ee ee eee ere me el per wee fh le byt keer | ere ke [ere erlergme yee) fe a en Oe CT ee EE Sere 
Be YT fies re auen Dae re | sae] a mel lalalalciaiel are ope peiene) er elem peers enmrpmper lo lere( | lel al | pelea (eral) | olotele (ol lape |e). (ol een ele petal Tt ty feet fete te Perry err rr 

H pte Best Nee dna eE RRRES EE RGGRnGeRSSoone eo A | ee tetas Slefe i elelete_ leleielel alee [ise ier i) oer a ele aie) ae falafel ert ei Grey tl ft ae) 

ipa |Top hae nolot ae pony alae eee eel Melee all) obser el ele me | el | sisicisteleiel | fale tt ii Tee ye | setae [| 
Fame Offa) Ff | |e ee oA] gk NN | a | ee Te Senn Fea cee | lene ts ats st | i pe ef et ge 


ae ee lela al Slates | feels alee ala aisiae | aeleiet ial tapas Tt tT ia 
(on ene ef ae ele ale alee ae alee ell ae Pele ee fe) 
elalsielele) Sejatelel [siaisialelsiafalele( fake se teetele laa lalate) ti 


rg] et ea ee Lape eae | ee ea 
eee fea es] ae i faa | ee alee] lola | le elle 
pelea eee ee lela Sls) Balepalaelea( | sale aa tae epeteleletaie fa) 
eles saielaiaislaelejafot ee elses laleaese tT ae eerie ie 
Pfeffer ol ae lien ata 
ioe conf pee eee ee ele ele eee ele eae eee aa 
Fa ale os] a Reape elapse aealeele eres eer Tae Eee aie 
ee See ee eeee eae eet ae a i 


SEE ee Cee eee eee ec eee eee eee 








Loy TT TP tT 


jelalels (eee Saisie CCE 


: eT esteem eet eer ttt 
ape | Sloe tel le eee ae ele patel late TT TT 








TA Tg PCA tia ttt tt 
a os ree ie er ee 


oe ae eae lola ie tae, lolol leo) ae) eerie) ae ee Pee eisai Ie eee ee alata elas Sa eee | de ey T | 

alee atoll geet Sl ae letey: pete) ap eae 

cee, oe ee Re elem gape ep lp aet Ile Rabe elo fo aera hn eRe taro mel 7 (mrenetal ale |mpmral vol ape) talele)a\a| slay. peta )elal alee eto alae remy! el eeeeete ll ier) jeg | ene) eee ea ela ee ae Oe ee eer eociatars weet TT 

fea eM ele [en ep ep ot Ores tele le) ome ela ay neem alae) aro) (ape lal alo lagmin| Spel Melaletal ane asa mpela alae) ojala fape ileal aa)s\s) an (see eee eee lejalala feel) a ( olsetaelliatats | Tear sae et Te Ra i ee eee meer ry re 

Be nb) oo eat ay a] Yea ONO fc ce fe g ea wet pn Fs TNF ny Yen ton] imac Nn | gs Naot caf cp memes oy en| ou gfe oneal amt fos] ony conf | peng] ee] eof el sel aN Lee ele re ale lolalmle fale) (ie paot aleeelepatale (lot fafelalalslatel Tae aie Peet Te aeel Te ee eee eer 
TPT TTT Tee Ty ee 





[cana a ee ee eal ee opal] oa elle else.) alee) ee aaa: | alee) s)eaiep ly a elate) 
LS eae | ES aire ean me MRM can ep ye ee fe ee pee a Ip ef lee) eee eae ae ol oe eee eee el ee aye lane el allele eal ae a ee 

ay te ae eee aie oa alU (menial) aaa yee, \elunlapel | | amatafe Siero foie | | Pepoto|: ajaleieio| (on eialafala pamela maaan) pee eer foe tel eer Tr 

fa Te ee es Tee lta ee Peneee we) eh i) ema) alee ee alae) [ate | alent (a) sal open) Peters CT TV ty | ee ee ET eli 


Seales see eels alee gee, el 

fe [ee ene has alee eres eel 

‘aee8 |) he oS fe eee fea tT | 
ga - ele ae) 2) 2) a eye ela sie) alae as es seiafelae ry 
es) Ape | ae ey TC ee ey 


ol ore ferret eee ee ae a 
See ee eee eS! ee enlen| ll ee eil lal 


HH 
a8 
Ph 
a 
Ar 
ah 
| 
: 
tj 


ple ere ele ee eee ae see eiaae el aeel la eee aa ee eee CCE EOE ie eer Ces reeieisisiee | i eer yt tr) ee ee sal ale el aaa e in 
ja ae] ele Paes ee aloteel sal pefelatey yale PT este ere Do ae late Tl het gee ee || ey lelslsfeseieeiet | See eee eer tT Tr rr ty 
ae a ee ee ee) ee re oa ep abatelok ®t | Ta] ieee eo) eieelee| seme eieat [ele Jel | letel siete alee lee eee ee eee ee ee eee: 
pea AE Ae © Bee eee Ope eee Leese PIS te te ee elas eke te epee Oe Bet eee eo eee neo) eee eee eee ee ne eee ee Cee ee eee eee eee ee ee teal an 
Tee eee eee ee ee eee ee ee ee eee eee ee ee ee ee a eee oe ee ee nee ee ee a ee eee eee ee eC al 
eo) oo el ele ei alalae sie saie eel TY 

pom ee ef ag eco mo] | ap 


ere eee Lh i a ZO 
Pg A et | WN Ree Ne 





pe th NE EP hte oH 4 — 
eat Cae A AN ARE Ee A Y  ed as e OCA A RAT OO eee eA ee ‘wis 
aE Leer i Teramae en pee cee) ol eerste eye ye iege | ere slepeet | alee (oi) lee COC Ge SEC SSC eee Eee 14 
eS apa Ce ere ee sme ole ey Tat Fogo |e Fs) ORO A mec | || | NA || | | ela pae) fale Nob tery fae tT ey 

. $s aon) oe afc] os fg om fms ff] fmf cada] aca J df 





. . L . ; jis) ane [ee fale yl ales elle lel ees Te yt le 
Lee ew te eee ee Ne ee or ee ee ee ee ee ae ee 
ot ce ee ae Br 5 ae ee ee a ee (A a ea rs re 
A J) es a a ee a) afm ea on ag pe fa fa fas fo ea hifi a) | peepee sf 





as om eae [cad eos ao yee nl Oi Ee Ce ed We ale me uel) | eee | Ley | ey 
COO TS ee ee ee ee 
rey el Pete ee ae EP ee 





Ct ee 7 Salta PE EE A He i eee Tt " 
Shetek as — = a —— — amet ae = a ar = 


—— ‘a ‘ = “~ _ ~ Sam a 





re at eo tts er PS 
Por 
Pt 


wi 
pera 








Pel opeapmpeee) layer g i fleet | tape 
SESS SEE ao Co pe ee rts 
a Coe ay ee | ofall fea a bet set ry ¥ Fm 7 nd . ns ae ms, 
PE OOOO yee tt ere +t 2 Serr Hor bh ee tT 
leo an ie) el we er ia ime) ole) fl leleeiaieien lel [ale elslelal seals) Ht toe bt 
lel julelet | lot we pan oro Tanne retail lat | mire mm | lee lame) lr) | ert ga pe Jer oe) | enel estetaiel bel bey | Cater Tal oe mua aes i) | to (a (eel meee intod [joke | foletlates | Ollie mi y TT rtry)ittts § 
ee Dp eT vit ih FS fH Lome A Oo eu 2 ee TT A per eg TT eet ote pees) eae 
FA amd LLY bed led Noel Lh be PIIA GV LIT Le ee | eo ev ree eae | felserey =|] =ia5 
ae ete eer eer yy Ht ++ Pe fosd [ala 


ee roe Nee eet ro oe eae oe leis Cee eee LE ee ee ae. 5 o aa Sr! = Ars. vet eee) lal 
glee)” Tele eer dense Pat TTT TTP crc LL LLLLLLLLLLLELLLLLLELLLLLLLLLL LE LIO-CL ET Ble S| | LO Was | CON VR Se PUA ae eeeeeeee—ee—ece—eeeeeeeeeeeeeeeeeeeeeeeeee 
iafmlmlnl wile at aif lal ale air tele lela lela elelake|=|alalaiel=|ein(aln[ets[o (iets) aac | ele taln [alee sieto] We) le) el_ atelier esto): (olor /6) 0.6) Se: ee eee ee ee ee eee Te Te OT eo 
aba nts ae alelieye le lie eee e lola lal Mel ate alma maak fol aimee alate ate)” (elaine |slaaetn lame, aleleleleter ye iaa|a\a a felaelmelee mere) emanate TT ert rrr rrr rrr Ce Oe Ee ee eee 
oleh SG Gie ia aan ian ieieieic) lelelee oe et le) ee eh lace! |” (elt (laren) ie | 0) eee) ie le ee ee ee ee eee eee ee Tne Tee 
PAT ee a ae ee ee Bi ee ee ee 8288/11. a We Eo 8 OF Oe (alae ase ee eee ee sla eee fala Sh SBa 

: etna ey 2 re | (ARAL | cbt ll Me tae Ve Oe ree ee 

et oe ey ee ee ee ee 

=< ei tt eeleaeeel aie el meee 

islapelel | lates 2) aes) Sal ele eee) la) 





Sp ia font onyan Lali Oa TC Fa fs atam [en ie Fonte at PON ses os fn] ah || ee eee] alo| |elalaleaial oelelelalalainpelala|e lan, = eleleleleleleret ale. oo) © ae (ane |) ele ,e tele, G06 Se .e la Tee ee ee eee et ee eee 
fore) a mare Se ee Semele Mee Lp lal le el yey el afeie | | fale tp lege Famke aff nt Ye] oO eee ee ee Tee ee ee Te ne ee ea 

Senior apo caret Fema | ee fae) tel alle mle Laem) ly pee | lela Exe] FFF USE Fon P NPs yes) ONT] me] wc Sem nF OME | cP eg |} fof ay on ned en een ge pO eT Lene py eee | ae Pe a ee ae ee Oe er ie ee ere eee 

CTT TTT TT ee ee TTT TT a] ota legate [lage tt aetey Sel alale| =| ieee seleer | ee ley olalaialeaie lo. Simla Cele elie) piielaeeta 


Sifmlarfaa hak eal tenet ale lente aT ee) a) pel ele eee lalate nmin in| aiela(eleln|ainlemieiala)sfalel wean | oc lelele |e) a peynls a). alee Ti eee ee ee ee ee ee Te ee re ee ee ee eee ee ee eee a ee ee ere 
TT oe ea) o) epee aaa alae. | alete laa lel elelaleal at eel alate tele. | le palete ¢ al” [ak 
cose See ee aie ee | el ae eee eles well aiee aes ola el fae) eta ey fe le aia) | 


1 eS RE Shee a eee oes) ele eee eee eee lel et Seale sea ealeel salah 
Sima ape ys rey ove ee ee al ool melo mmimias| (alee latelmiele lol m@ee|.\ Cael loop so. jogo) ale, Weletataeetalare| forelapel | Pletal wel po ae fatal feteletelats|agefeiet  ffel | Velelat y Talal | sieeieaer ice a lea ets feel eee ealels(s|[s(elels sale elas aaa eae ees eye) el saleae alee slay 
me Sees Lee | Dale eee NIN eee} A[M [Geb nels elmlepepes le [mgm lo) mpelmleiefele mama) | [alate pempalol. oy atera\ fo laler peyopepepepeape yt bepebe! Stole | set olen yt Pete Se steele TT yeisio TT eae ee. espa eiale ale slelelal aes a [2 olalal ala alalelalalelae alot (ol | lalelalaieles(alajaia( ial 
ts Sis sek Se Sak wel ae wee eves oh Beret ol | Foe lols leno lupmym epee ONC el) gyn) os) 1) [eae eee a ete aa eT lle ee ee lai is sae of) ra | ha] ey he | heey ea a, eae alae 
Penne ee ee cr eee fern | fee 2 eee ee nee ee eee en ee ee a ee ee ee ee ee 
Ce ee ee ee ee ee eee ee en ee eee eee eee ee Cee a eT ee ee “i 
ieee ee see seems (el eleyeyor fae) oleae ele alae merle fete aes (A eset elelataeie ie serie ee ie islet Peers 


Sais en eS Oe [lay ak Gk sare Jere) | eek a pees Talal lime) Et Soe tt ak ae Lee Te feeie) fa) pm 
Re ee eee ee eee eee eee ee ee ee eee nee ee edn nn ee ee Oe ee 

lel almielege €)epee ols (alae lel alase| selon |olsfaepa\ we la[e lene) alone(ala|e)afela(o)_\mmela| fale tatot (alo) Tn) Viale Aa ieee, Be Cl eeet IST 0 ne, fo) ete el eT eee ee ee 
ie oslo (ene ele] [aye] ele ahem Ole Oman lala [alalal malas elele|elain image (a). [olela(loje)|afo\n| ala ereyet test (alam) | oeeMe| oe) || ele age) |alalalolslaeleleliala(alaselsl |_| | cern naomi | (eet | en) CC ier Oe Se 





a Ry ee sn ke ee ee ee eee ofa ale oe oles arene eae eee lee siege, ee ee Sree ert lr (eis (cleteksletataletetet rt islet | Lae 
oh po + eleielatelaya(slereleatefelarere lah Pelstetele (at ate [weet atepeare erate ares) olerer) pote asa ee 


[eee ote) me |): ole, ee ee eee ee ee ree ne eee eee enn ne eee ee ee eae 

[afer atelale| stele) =| cme | Vote let at | sane ai Vee Tie ero lela o ters) elm—e a .e.e) eee ee eT eee ee heel bella 
A r {os eo] 4 oat cae yea Ge) RT NN eee ee ele, Se eats peepee ery  [olapeiay 

See ea eea eee lel see eee eee) eee ae eee ee iat) tater lal eal te) a alae ie ee eee TT Tr ea ae eee | Ciel peel See eee) pysleiee aes 
EEE EEE EE eg EE 
40 4 A PI s 4 

Paine on eee Been eee ea\ae ee alee e)n) a nee al alee salepeeeialame |e) | [oialoy: letelo( | oleae) ajale|a}elelalela)ojalafmleyo tate). | jalammeno_| Meiie eeal enim |e [ee ee ee 
= e : g lee | (stels|sats | teem et ei eisai ee eee 
Saelalalalaleielet taster Aet | ge) je woe TE ieee re 


fag | concn ee | NOD feng | fa) |e ee he tt 
(ony oy] ea ee | lala al oleae ee else ieee 
fae eee oleate lala er lereele oleae lalate ae 


COC ee eee oo ooo 
eA) See ae alae le eee ee) alee sees easel eee ea reel seep ie. (Tie) I oie iealers ater Tale Fela a eat ee ee aan os] ff a a) 

cee eae ae i tee ees Se lees sea ale. | [alae aele oieieleel eae eaapiet ey) Tan] tt pa et H ae BE epee eee) oa a ee lee [= a) a (le) pe cs jem aaron |p]  aenpete femm 
A Z [em eat] | ooo a | a sn ofa eles sll Re lelel pele ae eS 
se leone ae ee ae ee) sea es Sele Ses aaa eaieee) eel) lalaiel | eeeloaieai[) , : 

2] Set =) Bienes ee eee elec e aeolian ee eee (nial aerelalalatel | lalemet lier F RH et oer eerie 
se eee al) a a ee el Oa eee alam elelelele | ileal ib lelalec lola: eel laley ala 
eal): Cet es bake ale ae ee leo) 2) alee eae eee ee (ae elo ola ae Tay 








(elajaye] Je) Gaaea ee lol else ie ees Tee eel eee ny 
Fis) [my || || | || [| ca fap ao | | eels Pll Ameer lL St yee faye yt Te 








eine fi fs ase eee Tete fet Feeley Yep eae Tee] af fe eet de [a |e || ND | YT a amine miteal ola aleleletalanl jel alcalaletol=)| elf olalat| 1 fulona(][elapny) Joyal 
pa peso fe Te YJ Ji ei ff Ja) gpa ee | Dh | FO] fan a A [le Sela le lel alalapellale| ae, alolalepel || jana) | (ololatmiajo yay 
Fee ee aetna @eiiemigeco UL ioe Loe to to 





© ry) ae 4 

BE ee PEC 
ta ee po eee I ae ere eae at et Pe rye SoS gee 

a fy) eH Ee Ht + ay HE ED a Ro eae SD V4 EHR tet et 


[otal aes yal 
aT 


Ly my i] ies] ea ng nfo Gee) a) ap | | | fe) ee lee epeejeles| tT alate) eae sete) ape) eee cele ae eo aE aeitee) s/s palala eee a lalas| to ae, otal pelete, fal welfare palais 7 7] matte 
Peer eh eee Arh Stet ty Bete meee omg |g rpg ng we Yes am es | ae! alae see) ole shale) ) ley) Perel Te eo) Se olelelele a) ae alee lapels | tala lel isonet a joreielepelalope| | fare relay op) pa] : 
PT TT Tey LD lech lenpec dem Peed Tm en vinnie] Temps) ol as ORT YOY Tg ng ff oa) [Pf ss PNT ME Yea mtg || cp fg] ee Pea fen yenfes) Tepe. Tee Poets seleie le ele eleeeee | fee a eel sia leistetereree: (Cree lalate te tails loleeleto spe tlie) ey ere ty yy 
errr ee ee ; a Ae crt TTT TT ee A | Oe) Cs serie | ere) Sra a aale | laa ere alae eel een areola late iecweeregel | Flere leper. | : 
EEC SEER ESE EEC EEE ESI EEE CSCO eee ee 
Pd 
pt te itp el eee a rp ee 
[| PT eee a | ee Petree) | Lome toe. fe eee CT | d 5 Bo 
PECL Ee oe i CS Stott ee le Lh Talay erent ttt tt er ee Oe CE eee 
eee ee ie aa — Li Cte TT pap Ne ae ee Oe EEE UG Eee Tr ere 
Pre irteeLi L oat  r  ett  e  rre e bid) +} eee yee) Cy ter pare TY fale | eee ae fatale eT eel alelepale ls aietele | Taleo siete) alate a lalate Tie Cerra ately | 
Cle etalets| [olegepeyeleie— Jo) (olegelos e|olmie lege we [ol te| (oles | Viel ee eo ee core ee ee eae feta Sia errr ee et aaa ame ey Tare ae ee ea ole lala laaleler eae taleie te ee ete eS lero to fsa 
Poe leper slel am oloawe tate tome eet | Ta eee ee) err ey ee eee ee ee ieee aie) ae lee ae ele ee esis) 
ee a eee 
EEC EEE SEES EEE CEES ECS Deer ree tH See 
{ ss r 
a ah hn a FR ef 
~ Eee) \ | Q ad a PT eer ] 
icoeGeaeua et tels ro Reo ee ret Pott enh 2] A | lp OY a LY Lela tae cial) ees) ola aaa et ol alela = ele) cee eae esi a el Sarees ae) ee) | Tae als ae) atele ele laeapalaiaateteeteleymaiey sc fares, ret 1) farm amet Tela 
EEC ee ee CCS EE SECON EEE SE HSE ee 
seeeeeeees Feoen oseasoeeene a HH el) pate ates a afar) ng Yaa iE) eet sais [afealeefe) faa) japetae | ree) | ele ae oS oT B el) fZonceeea cia : 
fa) | MIN] gem 27g Rang Yen aN ng ge a Qos Let MM HOPE Ge een) nF] oa fem tes gape mae Pc a] ffi alee a] es eae ce eh [alfa ap seep eee eel awe cp aee eed: YS laalee eel perealay Tb ope) Me sale) |) emer Maialele |.) ame aie pelt folenn| | Nalanoqo ly sal 
[alee late eel apele ei toiaiepe) (epee) ate pee Tepe | ioe eee estate eae eee) 2) eaters) ae eee ae ee Saale ee ae elsiaper T Cais eeisia( ale stam fotatol | fueimpelo | ja 
a BREESE EEE EEE aR CEE CES 
= x ex i, ’ 
per epee ppb EE RR EE EE pen ee 
mt ia Be 
CEEOL HH eR PCC ee rest gener". rr] CL ee ae a a ee rt tri rier re Pree tert eee ee 
Fiale [onl alesis lapenel | falaln(al anaes enn Ole lsleae) (ol ilesf-| \u(ame) |otelopmlemml Gale! | Ng (0) fomeg sf ety TT pais ie jetereas ange TY [aie arr eae Oe tae lee ate se aes yan ee Saale aa Gr 
eet (Cae aie eae Belele () eananee Trt tii Lerriiri itty yey a eater) | tT Cr Tee eee 
_ plo oale sae este eee eee PLL eee lt et ee i en) | doen) ae) eet ry ee sees fate tok ae etal me fT fakate | oleate) flats | ee ale eae arya ay 
aloe) tal) of aolenyan one pele hoy) al ope mep cla) cial ee ceo t eae alae Spay ope h polato tee ogo samo. Posie ete ele eater elewiel | eles see ee le Tae eaaiee | sees lalsiape, [amare Yael) alee alelele lel) | epee oleimpern OP aEe | oleate amen o kat | | o] 
CECE CEES Co EEOC EEE a EO ee 
ple] ele) Lik | 
et taste cist enya ate ee iT ee sams eels |alalela ala f ate) alee alee ae leaele al lelalaleleia ele | late ale [ata a a) ae ela peed tt ele aes) ae eT alee) eiateiaeiatae) sale) lapale | oye pelea ar eee lett aint lal ital bel elon | fet) Ja) 


Le 
Fee Efe fa ea Pe fem sf fae ee reap peae YL sh een Ee] me pepe es ed | ne) Speier eae selec ae 
py lerelor=|_ | eiepaerel (etc) falake|o' aaeme | 1 (lee (fees (ok Pool (miaemie| Ponele tote! atta ee ee Ce See able leet ale ss slais sissy Bt ye ee alate paprveg ape] ey 2) Sepa pees eee apolar ee eee epee] ale | aero pmlael My sees ele e leat em) elara ala lele tea) 
ft a)mle| ole |afa) sla pe| elo lalole|a t=) eimnel lear fo eee) (oly peem mage ola fal mage RRL, [ge |ml pea fe Es) olen fepm peep mye pe ey TPS fep pay Tf lames apa Please) ee ape ea eee) ae lee ere ree ee eel a eels las eas aaa) Salers eae ese) e) age) ape slo le be es ae ofa lapelere sla pae me gate to) la 
by enten Of Sale) meet | nee cele etme ft eam (Oy mmm mal egne ie lepers lalol la) alae aaa aleeta saya lol slate) pee elelele ae leleels aateeelelsleie| "(sae le (sale e)s ela arse aaa alee) elalel als ola) ease lal slelels isa ate) altel tates! olelefetel latest | isle a aele eee e|ebs)alat 
pf ayelol loyore | epee Coe ey] taper tee ae a Trt eC CL told bolt ne et) i al eee ot ime) i le me elt) | po () (open Te fet te) ot tet | | hr) tet yo Pt tata) iene ghee cio | tote 
Pf ae alee lapce mee ap oletal | el | | ee ele eee Sete tenes | fata, | oleelale| laleleieletsielepe) lalate (pales) leleleiatale sale) alee lean leila elem taalel Te eee) Fale ake alge fib at] ee ole yey etal Armee Cea et Ala fone fehl) 


rf 
— 


Tar lett cla(aneelapael | ale] | [oman “year oe aoe S[Ole) | f stemmiel | nl mpmpalagete | @errei a] | 1 |) alepelelaia)ajael am ympelag ake) papay fenel | Pumper aeatefmy a) ajm|e pea m| amie) el lela e lear ene elale| feteeeels el alele) | aaiemmele) | | lame) (ae ere) lel maln) foetal! [ol el fermion © Letemteter@l | | amet) | ojape) s(n appa amtain(=}=\e1e) 
pal (logo) | slepapmrelaie| | a(a(epalal pel: sal alate beef | feiaelelalalapa alae) uela, pata | ata alae eee apa) of ane) | eps) ola see aleele yells allel oleae ele elo eee ela ey areola aloe ae as ae elalele ol lope Palaie pera yatet a eee le pale pak laa ays late f a) olay Wie a eee | wpe 
fay eee iaelete! ae ae go elas ae ee eee ee ee eee eee ee ee eee ee ee ee eee eee eee ee ee ee ee eee ee ee a ate) @ ele Geel eae a ety) 

[atoll alates fey fee [apne ga amin =| alarm) aol ae Cale|el ml eim aga) s ae ee len (ale) a) sale|erm af alerele lela le|m—fo|a|elafa( (alana) |apmpe)aga\e myo) [olay mal alates) (ane) alma (amma |e | Nagel” | fo lalula|m)o fake mel | o)m gm) ca) nla] ale] |a pm imia| aa pm in! 2) al etenl aloe lmlala)s(slora)-\alaipia(a|e(alelal= |= (at fo) _japape paler) Nelo po ote ya mela 
PRE el aes ef ae alelele alee) leleteie ss atelele) io) aie see lee lela) pt | aS EL pi ppm Pe ee | hale se re Tales or ee ere ene aa a 


foxy pai aap a ae ah] 

ee) aye epee alelee eee el Seas alee ae S| laa) eee eae fale a slere)a)apa(al alae aele, ae) alaae eee! elapeee lope aes) To elatelon elaiepeieelae| Tei ta a feel) etal 

TEC Terr a ee ee ee ee ee fate eels) fale lela sale lel alars ale laleeieletel) | le) 2 tee alee lege | oleae eee a) eee eee eal aiae eee ease aeepelapemet | Foletel | Petia stat ieee en alee etal 

PTT rrr rT a CT iene TT rr rT eae a re eee aie fae lee eel | Pale ele alae ae elmer ee) tel el a oes) betel 7 at aie Cyt is) 

Plcmcloie lope | lorspol_| | hs | | Come |g ate ee pa le peiapeieepot | ee Celle fate wa [ly alt |) Talelelas elle sieiel | lalalels lel ae) ales aerate lame Tl tele eee Flaite lala late aateee sielelr [operat frag | yf ater feketestey | late a a Be ey 

sqaeereat | Vlopway lero TT 0s pepe | See ee eee allele t lates peal te eee alae See ea ae eee ee Siete ae lee lalate Tegel see eae lel fe lelele alee sl olaieleeniate eile sles ee lies | ae ee eae eeieiser hy 

Be ee ee ee ee eC ae ae eee eae ae Pf elope Pee alee le eee eer tae are ere | Delete tere crera 
ITIL Tat iA IIT rr a ee eee et ele ae eee ier alae) ° 


ee De tec impale le ee ele eleieaie) las ale a eaten ale ote ae eae ee | eee eee ele eer al leet ae ere ale lela el ete et alee ete iaiafel ret ieee Ee alee ister y Fo 
PTET TT rere ee a ae ae ee ee ett tt ayy Teta | alas eer eres at arate a tele lsiae| eee) ae lelel alesis ae ae ee alate slate lela ppt rafal ey seme [S| ape | eee eet wl) Totals ls ialat ll alenme pete y | yy 
Fahy ia) lat ee eae lela aera Cael) alate alee eater eee ee ae ee ae sige sleet ee ee Cae eer Oa ees Ie eee iat ee ae ee ee a 
AO OC ae ee ee ee een) en CL Cee ea pe ee) Oe ee eee) oa isi aL ae Te EO eee a Tt Deeaseaeee bedi A delete Leia altel Loa 
By 1) lel io ates (i et Ce ee a ee ee ee ee peta rebate gee “a 
ee ee a eee ppp pag See) ane [olan eee ee eee ren er ial lt eee ee oe Celene fer el Tr 
rh me ef fem ang avenge pg Gp heey Vemma | peal pam esp fs OPM on) og maa Ng gf no og ge Fog 8) fn] fa oa) co Kee lg eye a] SP a pCO YF ceges ey ep waa on Yea ap) falafel aaa fesse Taal lelpel al elenee) | alee Cal else ie ) alae ssa | tale | baler Lee ee iat Tre tar Ce 
ey) ela pee (alate o fala (ale), eee | meee thoy | folamey | foie!) olete eo iaiel alesse fae saleteter Ts) Fae ee ae CS ee eee Ei SSeS e eit ere a eee ae es Toe) tente) Cal ate Tae rT See re to 
ty TY oiataleielsislelapel oi [oe | eee | a oe ieieie eiaietse(atatsolelejetslat | late letse | eee Te etefeleel atelejela| ee efeleteiee fel teleesietesle)atene ete la ee ales sails a ae ae) Se ee 
PCT ee ae ee eer ee 
Piast fer islelelaleatebe( | ial: [eee | eee eee] | Laie eee tae eleleiatelelelelatapeiets) (elalaeieeiere | ee. Biel feleis slate [lelelaleyeleg el shepetals ele aielele lalate lela laiape | laa fale eae eae eS ae eae ae ee eee ae eee Ee ee eet 
POL CEES Ea a gelotelet | Cie) ee ee CEE OCG OE OE OE Te EEE eee Se Ce eee a eee Let Ta) fapey pea 
Cre era ZA ae ‘eqs ECS Ce ERS OO EE ee ite aise sels) seca iels) eo et eel are ae eet Tae le ets eee ial al Tt a Lich Bethe ef} tt 
De ae ee aes HGS GH EE Ree ee 
Oe ee ee ee A ee ee eee 
Bama Lee fe fee og fc Jara pee | |e) | eee ee eee eee ene ol oS) eee es) sala eee aaa ale lal sels alana.) ae ee a ee ee ee eee | eee ee eT ee eerie alee ae ee ee a eteere lala 
Bi a) Feeley aba gg |] | eee ee eee) ees) (ote) et | et eae aeiae Ss isiaiet | stele lel | sate]. eiee cla islsla(afalalel ae) seieiaa siete Ce ee ae Tae Sel arta aa ate | Taye] by pp ae al els alam eee ler) alee Teepe Te 
pe oe ee eles alee [al aaa ee Ie eee) Ge ee a a ee EEE eee eS ee a aera eae et aa eth gel ae) ae aa a a 
Deel os alee eos alee ia ieioleee | tie) | el eleieliaeee | ialet Maes eee ee leleiaetey aise Te eC Le aa a a es efetateaiael_ [see aie) aie(oy ae pat eT al ot eS eee pe ae eee ae 
ee eet ee ee tee ee ee 
Pee yells lag te] teks) slaee(e(elel=)alere eae alate, eee te la ee ie) | eae) sf ae ese ay te Rees la de ep ea es fee) a ee er ete 
ste et eee ieee eae ale ae slaiete siege sete) Coie TT el Siete) 1 eae cl eee ee) Patel t sielel fapeare! | leet eer ere ae ee las 
ea fae as) ee eee ce aoe rT | palate CTT ee ee) eee) a ee ae eT a ee Jot | ae 
eee eee ea oa oe eee eee eee ei hee eee ed ee eee ee ee re ee 
Cee ee ane eee ee eee ee ee eee ere 


OE STDs fe ee i im ga ef ae a ma es ema a fs fe) eee) ae fe] pee es) fee ee) ee eal eee ele erate fae eee ter) | alee elo. ee ae S| aieeee lee [ema (aia ae a) 
a Sh Feel ea Be Tai fe faa ea fea fo fea fs fafa eee aha eles) pare eee patel | tee eat eee ete | ee ae oe is tails slelet | ieee fletels siete sitet eelslalsleials) eee aaa oer age | Laie, (isabel ieee ee eae ei ala 
ear Ceara asf Tea ee a aca ap eg age | ga ee a eee fas) eee) | Poet | ae) eel alee te eae a ea ee Gelcieaaeeel eelels)alaee) Telefe Tape le lapalel alae Tanta elee! feletare | ase ae ee ere ee el 
se een) aal eee) sie alae |S) ale alee a eee ee lee lee eee ae | lee) etal aiaalel | eee aaa eelaee| eal ial 7 See eee ala ae alee aleelelelele ee sl isiel ialsielaleh alert | lo faaieaelera tale) Vata aay Sete eta elt ee ae ee ae eta al Sl 
Cs ai tat le lanelets alate eae eee (aisle le) a eee fale) olelel aie eee Ol aleleiao alelelelelarer (oleyane) | a cal o_o eae ae eee ee ee ee ee eee ee ee Te ee a 
ETE Tee fee ea is Fama ape ee ee | aD se Pam fs fa | fe fe seam |e) ee) ole falta al tale spel ere weet | el ell ralel rere) |) | lame | ate) To lelomoe tae ears) tole la es eee age aise leleee) ae as el ae ee oe a ese) ae ee ee te 
pee ae a aee eee [alae ele sala lelelaae lela) alae alee eeee els alae alae) ss paa eay eles (alele) ola alee sale a eee eee ee ee eee ee ee ee eee a ee Cee ee EE ae eee gee reer rs 
ete ala oleae alae epee) ale lane) | lelelelee |= eeelelolsale| aa epelalalalo me eee alae al ere ape Ty eter petareret | late, Ty ae far ae te epee es 
pea ale | leer seret: | late ts) tae le ela eee ee et 


Cee ete r als) fe Seta) ae ol PT TT PTT ad Bh) 

@ @ j 

3RSRRas 3) i aya ai | folie Lo sla Te a as) aa |) eae eel ian | lar ae) ete) ie hella] | la leary ota. | T_) ae eels) oo lele(a)r ieee) ee alee ake) Cal 60 ee ee 

ete) ty tl aa sale) ore) halal L Bo - F ai Pia) oto) alae) Yl att Toa is I T To) ese) | ea to) flo) fal ae) |) nel a See) ecle| \aley Jo) | eee) a lor | et Pe [aa | elojae TT pap) 
jit Ty bi ae! s i. a8 : i tL 7 7 alla oll as aa) oe) [ee || toe | inom | | moe.) ) ao | 5 If be £Eese 

gu - P ret Fr 7a) aaa Papen ee) , 
- F ¥ “1 Lt “Ett TalyT rT iT Ly od 
ooo. Tis Tee eee 





abel ila Tas teeta etet tat Bele) a) alae eee) oa) alae) aerelele) tele 


Ee 


— 





: fl eT Ta my ane TT T 
a HE ue a 
7 ea 


Vee 
a a ae 5 
LL _ 


al aaa ae i Ce SoS gay a 
| Teele 2 : | Le sv Aa Es 
HH tH ne ae HT [ ih PSI ak if A pf HH HH HE ie a ne SHEE HH ee ah 
ay 


: He HE Pace Et e a ti AE a EE a ae iE Te FLEE Ht Hie a i HLA 
ae al A ET 
ae a HEE oe 4H FEE 
ane TT TH UL EEE HTH i ae eee 
Bd a | i cc HE A i a a LEE en 
ee 


a 


ESE EERE HERR LE TEE HSE 
He HH a Te I 





EHH BEE EE 
VN} a Eee 


t a ne ie ne ic : i ede a iF ae Hc ( ee aH CRT 
mae ti a HE Hate a iC a Le VATE He eae iat gts SEE 
Ba) HH Hn HEELS EEEE EEE oo 
aa if HAE eH ey | HE HE Tae EEE HL: FEEEEEEE EERE itt 
We He rH cree FH de cole PH Ee i He ag HE - ii he | 
HEE || Heed i a 4 a Lif Hit EEE HE HL HEE HUTT es 
ae HTH \ He a Te 
a 
SHE i ALE He 


AE ae HEE a ai 1 
THEN oie HH i i 
_ | - Ht aL aaa ae i SEHD a 





PEE epeeee cree 
oo pa ee ee = 


| if N\ wa He a HE | | 
Beet at 7 EN Vi. ESE HH ne HE 
PEE et ee HEE HE ao ane HE a it Hat ae Me ae HEE Ht 

NEE PEPE EEE eee 
a CC 
Ce TT dae St ; aN VR el fete ae BE 
EH LEE Ha THEE vue HALE i A HoH FS Tie SHEE eH 








eee 
. a GT ee A ae Al a a rapper re eae stat i [i 
a ian TE Li ae Ae ne ae ee Te 

BUH or ae SIE EEE 
i HE ie aN is . HTH, i HAT iP HI HH A tt wily Aili re Hy Hie HL i a ue 
a fe a | al i 7 _ a i Heh 
E a an ‘ ATH AEH Vf TEBE HELE Hf Ht ti HEH HEH HL HTT neh HHT Ht 


ag cae au | Lt i of / te a : a a ile 


ei 


HAL, Hei HH fi EE ie ae dE pits EH EHH ie aE HAVRE 
PE HEE EE HEEL HEE HEE fil HEE HH HEHE HH HEH ae oo 
Se i 7 eA a eas it He RE 
mee ae FT WHEE 


RATHI THEE HEH EEE Ha FY He HH ait Ht te 2 Ht He ath sit ri Hi tp Ha ER 


a Ha le in - i ae ul ci i Te ee oo 




















el aie 
ha | . Meal i HL EE aa tet SARA Ha LT 
EH iit HEE Ha Fp ci tb vate ay ane ase ee 
TE aS Lite et oY per ETT Cn 

Ht cal FE ate " He HHS EE ia ara i iat Fs 

ie ae ie aa il 

BEE a 4 

PEE THEE PREECE AHL HALE EEE 


ee 












































ECE ETE The Hr 
HEE HE HH A Ht HEE PEELE HEE 


A 
BEL ae AL TEE THAI VE VHEE HIE ATE HEH SEH Setieatuites 

Ha HEE He itt EEE TTT THEE ANNE HVT Ge Tee EEE HEE Hil EaHad HH EEE HR 
| FEE TSU Ce EE CEPR EE ener 
it i HE ae AREAS THEN Hi VM SHE HELE aan ESTEE 
oe FAUT EEN ENTE Hie Hea oe a ar EE ae 
Ts ea i 
a Heh bee Ha EE He HEE ee ve a ENE NE ao _ a ie 
















































































an EEE 
BEE aa fa tia ial fee Het a 
SEE TT Tee oT ee Tee a eee ee a 
Berens UHH RTE Pee if CHE Ee ae He ieee Hi 
leet i Cn RSE TT gat eee 


























Prete 








7 \ i 4 i ee 1 { +4144 en I 
Perr peer ao PT | 
ae dash Tplial Sua c fa met as S| og | 
ge ami en Ye Fa ag a | ee bn ad 
SS Fea Fee ea as fa ee fens Fem ws Ff gee es |e Fe | sg ff Fem Fy Oye fof oot | fife 
Tafa a) eae a a Fee ee a Fe fais of Fen fe Ope ne FLT 
fanaa ge] ee eae ace | 








| , we 
H 




















Hap py 








ine 


PERE EE CEE Herero ee Ree eee ree rere eee eee 





i 
Hy 


EL 
lt 


pL] a 
a ae |) pee) | eel eae potas ! TI 
HH sans He VJ ere poe ey foe eae) SouaeseceeenEans se a ae aa prep 
azeaBs egies) fe] ha) ete ose ee see esl ale lees) mm Ff a gp jee ei Teele ae ee aor [Se iea ee ee 
etait ee aaa eee a H H eo fa Se eae a ees a le sey Ryle Ta BESSSASIRGPOOCOITeA 
ele ese aa [a es peat | as 
hme fe te | eae | ee faa ae tain lml as 











pats) Sle (alles are) a a ey fam i Po se) eee | ye 

ee ae ee eee el a eee nn ay apmiale, (alaea(s) etme) (laimetemt fo) | jatot | jateenn| 

Ceol fee Se | en gee ae i eee lele lle tll eee 
Bea oS eo a fo genoa a enon ee weg (eT [| f= [can] PT | a | 


TLE 
I 
i 


it 
HA 





pat 
oat 
LEFEH 








fH 
) 
a 





BEES CEC CEE ee 
Seen Ae ee lee ot ne) en en te ee TT tT jot T ae oe) aie aa ert eee Biol | ee ia ee ae I leer 
See ee kt ole keel lara pte Se eras | [peels pe el ye) saan a eee alee eeepc m lol lage | ae fellate |e fatale) Lemme pes BURN Ss Lem efe| ol ef opea | epmee fale lol (a) aie] alee tele ves ere open ee alor) lan paren a 
BEES EEE CECE Ee eR EERE CEE EEC eo Ce 
ee oot oe ee iol ape tet wl | eI T ti bt) Cerri tee : 








TE 





ae ee ea ee el eae eee ae 
eT aU eS 22 a YL ita ew EY Eee cae fn o 
i) ee Cele ate lo olen 1 rt |_| ero ee, ete 2) 1) elm i 
Coenen ee ee eae ae er sp Poe ee lone a) ae ee TT tt lee yt 
eye ae) eee aaa eee eee Ty 
Cae) eee ere | 








jaan] reas) Polenta eee ea) rere 
JE DDPaeaen 7 Salm aimlalalotalaletets ictal aes eae) sae eel ee ae el fate ye eo Weal vey fy yale) ye) stad) _) a 
OTT tt jt pt — th eh alee) +4 tt te JOBRBP ase) ae aig) erslel aaa tsa ales) ar et | loo opel | Joelle!) ite l=) ele ae) eal eleanel. (| laiee) | 
booth +t Soe i4- Sean eeeee ee erate aaa) 2) alee) els fe)e lake) ale aie a elaleeete) ee SOC aeieae tt 
a L. J i I SRM ae Ue: te ae ee ee ee Lo nlet tL tenet? | tl) Toei) Sate we dm fa nd i th 
ait T fics ais + 7 Lf! a T 


4 I J ie 
J | é 1 i ‘ - 


eet eee a eee a 
A ake ee ee a ee) ae 








Se | 
eee PEER 
- ON ce Hefeeseeerte He BRE EE SEES ee Ee EE a etree 
Peer Peer ARE “hap 
NT HET AEE EEE 


BENGE Eee SEE EEE 4 ae 
FECES PEEPS eee ECU ere 


ee et | Te 
i 
i 


HH 
PCLT 
Set 


lo 
zi i= 

i 
an 


| Teel ere 
cee il 
Yi | aay || i 
{ 
{T) | 
ri tH 
Wp 
' | 
[ai 
{ 1 | 


eo 
b: 
\} 
ri 
. 
1 ; 


SE 
me 


q. 
,—- F 
a 


TER ee 


eI 
1 
ee 
satay | tt 
1 oa bi 
| 














ae a 
: H+ AEE eS fe ee JENNER BEE 








: OC Hoe 

a ECE 
: 

eer 

OEE 


yim 
amin (2} pe bfe| 

¢ i 
La | a 
I ' Vda! 


ia 
TPES 
. f 
i 

i ae ee iy 
fey br 
! a i 7 
rj ; f 

{ ra 


tert 
a Aa 


























oo 
ae ae HE Sea 


SH 
FE 


\ 


fa) St ae 
gt 
ik 


| 
£ - 
i 

















mee ee 


a 





ETT 


pel 
nn 

r 

ae a 
ATUEUILTH SNR EL als 


fay 
—- 
® 
& 
a 


td 
ah 
Pt 
ie 





nun 
He 
HT 


nt 


i 
A000 


3 eee Se 
| 


eH ae 


+ 
eat | coo | Li 
ee aa ty ; 


SEE 





| CPC E 


i 








C. Evaluation of Probe Measurements 

In a dense plasma, in which the mean free path is of the 
order of magnituce of or less than the probe dimensions, strong 
deviations from Langnuir's theory can be expected. These deviations 
depend on the size and character of the probe. By using the 
approximation descrited in this thesis, one can caiculate these effects 
in terms of an "effective probe length" and e "potential disturbance." 
These two parameters can then be usec to analyze or predict the 
experimental results in the manner just employed for the proiate and 


cylindrical probes. 
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PoGnetiCetrhel> END EFFECTS 


Let us now investigate the effects on the particle density 
distribution functions of the decrease in magnetic field intensity 
at the end of e coil with attencent intersection of field lines vith 
the contginer walls. The effect expected is that diffusion aiong the 
fielé lines, being greater than that perpendicular to the walls (radial 
diffusion), wuld cause the density distribution near the ends of the 
field coi] to decrzgase more with increasing radial distances from the 
axis than precicted from radial diffusion alone. Particle losses by 
longitudinal diffusion to the end of the field coil, where some of 
the field lines intersect the container walls, would be responsible 
for this effect, which would manifest itself physically as a con= 
Striction near the ends of the field coil. 

An exact solution of the continuity equation for this conficuration 
(Figure 5=l-=note the enlarged radial scale) could not be found, anc 
consequently the approximation described in Section IT]=A=2 was used. 
The plasma is civided into zones in which the diffusion would be 
primarily radial or primarily longitudinal (along field lines). 

The latter zone, in the vicinity of the intersection of the fielc 
lines with the plasma container, iS 2 zone in which the particle dis- 
tribution is controlled by the loss of particles diffusing along the 
field lines and recombining at the region of the wall where the lines 
intersect the wall. The size of the region of parallel diffusion 
increases, of course, with increasing magnetic field. If we neglect 
the curvature of the fielc lines and use the boundary condition 

aoe (2=99) 


ent 


102 











we obtuin for the density distribution along a field line 


N(Z) = 7’ sun I&, 2 ) = Wen(%2) (S=1) 


where n* is a normalizing factor to be determine?, the z-axis is 
the cylinder and magnetic field axis, and the origin is the inter- 
section of the field line and the wall. This solution cannot be 
rigorously derived, but is justified as an approximation consistent 
with the accuracy cf this calculation. Comparison with the solution 
for a bounded cylinder in a homogeneous field, equation (3-62), shovs 
the same z=-depencence for the density distribution. 

A second zone in this plasma is that in the vicinity of the center 
of the field coils in which wall losses due to radial diffusion to 


the walls control the density distribution. Here 


~ Fee ) cs Mates) 
FOr y= 5 i Dm) = No nt = im (S=2) 
where ny, is the density in the center of the cylinder. 

Then, as a first approximation, we shall aSsume that the ionizae 
tion coefficient, y » in the zone of parelleil ciffusion will not ciffer 
sicnificantiy from that of the corresponding Stable radiai configuras 
tion. In higher approximations we would have to calculate the change 
in this quantity due to induced changes ir electron temperature by 


the different fields, but we shall assume these to be regligible. 


Consequently , 
yey (S=3 
Therefore, ve obtain from equations (5=1) and (5=2) 
zy. [Dm <7 ( oa 
- Day IN 











The surface separating the two models (equation (3-l1)) is that at 
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Substitution from equation (5-l:) reduces this to 
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Figure (5=2) is a nomocram from which the zecoordinate of the 
separation surface, Zz, cam be obtained as a fraction of Z9 using 


r (as a fraction of cylinder radius R) and 


Dm, mn 
A = [Pe = Yl+ pp (5-7) 
ab 


as arguments, From equation (5-7), it is apparent that “a” decreases 

(with B) along the cylinder axis (indicated on Figure (5=1)) beyond 

the end of the field coils. Consequently, when using the nemogran 

to determine z,, it is necessary to use the value of "a" at the 

junction of the two models. This requires a trialeand-error precedure. 
In order to illustrate this method of approximation, and to 


estimate these effects, a system with the following dimensions was 


selected: 
Plasma container radius: 205 CMe 
magnetic coils: width = 15 cm. 


spacing = 5 em. 

number - 10 

inside diameter = 25 cme 
outside Giameter =- 50 cm. 


overall Llenath = 195 cm. 
10h, 








| 





The axiai field of the above array of field coils vas computed with 
a dicital computer (figure 5-3), using the thick coil equations from 
reference (20) . The author states that they agree to within a 
few percent of measured values, and for the region within the coil 
array, experimental curves from a very Similar array were available 
and agreed very wll. The fielc was assumed constant through any 
crossesection of the cylinder, which is a good approximation because 
of the small radius of the cylinder. Figure 5=1 shows the field lines 
as computed from the field strength data. 

Figure (5-1) also shows the separation of the tro zones of 
diffusion for three different orders of macnitude of ae nanely , 
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Be = 10,000 -! 


spite fl - ) Ch = 100 
(5-5) 
It is to be noted that for B = aa (a =10) the effect of the 
parallel diffusion is very small, and is centered at the end of 
the coil array. For Bid. (a = 31.6) the effect is appreciable, 


and extends more than halfeway into the cylinéger. Figure (S-l,) 
is a density profile of the resulting particle distributions Figure 
(#5), selected cross-sections from Figure (Sh): and, Figure (5-6) 
is a plot of constant density curves, 

One feature of the density distribution rhich must be discussed 
is the question of determining the density distribution of the central 
zone when it does not extend to the wails but rather terminates at 


the edge of a longitudinal zone. A longitudinal current (between 
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electrodes at either end of the cylinder) is responsible for the 
ionization within the plasma, and this current must be constant 
through the plasma. Further, we have assumed that the ionization 
coefficient, D , is constant. From this it follows that the electron 
temperature is constant, and therefore the longitudinal electric field 
must be approximately constant. But if the field and current are 
approximately constant, then the total number of particles ina 
cross-section must be constant. Further, the plasma in the radiale 
diffusion region vhich is not controlled by the cylinder boundaries 
because of the intervening region of parallel diffusion is expected 

to maintain its Bessel distribution (thus satisfying the diffusion 
equation and not requiring a different ionization coefficient). As 

a consequence of the constant field ani current restrictions, however, 


the tctal cross-section particle density (particles per unit length) 


K 
i nents) Me (5-9) 
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is a constant, equal to the total number of particles per unit Length 
at the center of the field coil. Therefore the normalization of the 
Bessel distribution for the center of the plasma in regions where the 
distribution at greater radii is controlled by longitudinal diffusion 
must be so determined that equation (5=9) is satisfied. 

One effect which ve ignore in the previous formulation is that 
due to the small radial component of the magnetic field at the end of the 
coil. This component is strongest at greater radii, and decreases 
the mobility of the particles in the longitudinal direction. Con= 


sequently, it would contribute further to the constriction of the 
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plasma at this point. We neglect this influence, assuming the 
curvature of the field lines to be very small (recall that the radial 
scale of Figure (5-1) is greatly enlarged). 

It should be noted that there is a discontinuity in the directim 
of V n along all the dividing surfaces between zones, and a discone= 
tinuity in particle density, n, along the a = 31.6 boundary from 90 cm 
to 130 cm and the a = 100 boundary from the center to 16) cm. The 
former discontinuity is a consequence of our approximation. The 
density discontinuity occurs because the dividing surface lies along 
one of the field lines, and is a result of the fact that the sine dise- 
tribution for the longitudinal model cannot be made to match the 
normalized radial distribution along such a line except at the end 
which 1S nearer the end of the coil. The transition regions are 
centered along these divicing surfaces, and are the region of greatest 
uncertainty. Along those surfaces at which the particle density, n, 
is discontinuous, this approximation is even more uncertain. Because 
of this, the a = 100 description, with a discontinuity in both 
and n along its entire length, must be much too crude to be valid, 
and further attempts to fit this approximation to the effects of 
such large fields are not attempted. The description for a = 31.6 
is believed to be a valid "“order-ofemagnitude" approximation. 

In Figure (5~7) the heat liberated at the walls by the ions and 
electrons carried to them is indicated as a function of position along 
the cylinder. There are three curves; all calculated for a = 31.6: 
(1) the wall heat loss for normal ambipolar diffusion without a 
magnetic field, (2) the all heat loss for radial diffusion with a 
magnetic field, accounting for the increase in radial diffusion which 


accompanies the decrease in magnetic field outside the coils, and 
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(3) the estimated heat loss for the magnetic field, taking into account 
the longitudinal as wel] as racial diffusicn. This latter curve is a 
very crude approximation based on the resuits predicted by this model 
with adjustments to compensate for the errors introduced by the 
discontinuities at the separation surfaces. The figure shows that 

the heating at the wall within the region of parallel diffusion and 
within the coil array is greatly decreased because the particles are 
transported to the end of the cylinder instead of directly to the walls. 
However, outside the coil, where particles from within the coil are 
carried by the longitudinal diffusion, the effects of such additional 
heating is masked by the heating due to the increased raidal current. 
which flows because of the lower magnetic field intensity. 

The values of magnetic field B for which this effect can be 
approximately described by this model are those whose order of magni= 
tude corresponds to a 431.6. Smaller fields (a 10) will not cause 
an appreciable end-effect constriction. Larger fields (A = 100) 
are not describable by this model. For a plasma which fits the sample 
data given in section 3=-Bel, this is seen from equation (5<8) to be 


approximately 


QGGOY j 
Ba /Lt a tb weberfm2 (5-10) 

7G 
Tt must be recailed, however, that the above description is based 
upon the walls being perfect insulators. Electron diffusion, pre= 
dominately parallel to the field lines because of the small electron 
Larmor radius, is restrained by the build-up of a negative wall charce 
in the vicinity of the magnetic field line intersections with the 
cylincer walls. Similarly, the ion diffusion, which would be predominantiy 
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rndial because of the larger Larmor radii, would be opposed by a large 
positive wall charge within the coil. In this manner, the conditions 
of equal electron and ion radial current ami equal longitudinal current 
are met. However, it is known that insulators have a surface conduce 
tivity which depends on the material and the temperature. The effect 
of this conductivity is to permit an exchanae of carriers along 
the wall. Consequently, we obtain diffusion which is somewhere between 
Allis diffusion [15| and our quasi-ambipolar diffusion, The wall 
charge buildeup is avoided by the exchange of ions and electrons 
between negatively and positively charged zones. This exchange results 
in reconbination along the lenath of the cylinder and would probably 
invalidate the conclusions of Fioure (5=7) if the currents were 
appreciable. Further, if this surface conductivity were sufficiently 
high, the diffusion would approach Allis diffusion, under which condi-= 
tion the ratio of the parallel and perpendicular components of the 
diffusion tensor would no longer be 1 + Meee. but would be 
1+ 112 Bx Thus the field strength which would produce a con- 
striction describable by this model would be reduced to 

Bs page © 0.4 unboe/mn* (5-11) 
The description of the density distributions would not be affected 
by the change in diffusion process since the only effect on the 
diffusion is a change in the magnitudes of the components of the 
diffusion tensor. In practice the diffusion would be somewhere be= 
tween the extremes of no surface conductivity and perfect conductivity. 
Therefore, the value of B to which these calculations correspond would 
intermediate between 0.) weber /m@ and lh weber/m¢, depending upon 


the surface current conducted by the plasma container walls. 
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It should be noted that this description predicts negligible end effects 
for the plasma described for fields below about 0.1 weber /n®, For 
particles with higher mobilities, however, this effect could occur in 


smaller magnetic fields. 
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APPENDIX I - Curvilinear Coordinates: Transformation Equations 


Because of notation differences between references, the trans- 
formation equations for all coordinate systems employed in this thesis 


are presented below’ 


1. Cylindrical Coordinates (r,g,z) - Figure I=1 


x=r cosg 
y =r sing (iam) 
Z= Z2 

where r&=0O, O< 9S 27K -00 <2Z~< 0 


2, Spherical Coordinates (r,0,%) - Figure I~2 
x = r sing cosg 
y =r sir® cosg (1-2) 
Z = r cos 


where r >O O=9s527 o£ 06T 


3. Prolate ne Coordinates (7; 23 d) ~ Figure I=3 


2G -y NG |" “Ce 


a Ly Yes yy’ “sD (1-3) 
Zs 3 1% 
where <1 <1 < 1 1 = <— © OsP7 < an 


In the prolate spheroidal system the surface & = constant >i is an 
elongated ellipsoid of revolution with major axis of length af and 
minor axis of length d CEE The degenerate surface ¥ = ] is the 
straight line along the zeaxis from z = -3d to z = 43d. The surface 


71 = constant <1 is a hyperboloid of revolution of two sheets with an 


ine | 





asymptotic cone whose generating line passes through the origin and 


CUM ake 
1 


is inclined at the, (2 = cos /| to the Z-axis. The degenerate surface 


| = ] is that part of the Z-axis for which | z| < $d. The surface 
% = constant is a plane through the z-axis forming the angle 9 with the 


x,Z-plane. 


h. Oblate Spheroidal Coordinates (y, es g) - Figure I= 


x= AV eel” coop 
* alo) Ug+ Nl? sng (t-l) 
Z = S43 


vhere -1 5 7 = 1 0L2@ em O<f7 £27 


In the oblate spheroidal system the surface Ie | = constant >0 is a 


flattened ellipsoid of revolution with major axis of length aye Le 
and minor axis of length | 3| « The surface 2 =0 is a circular disk 
of radius a = 4d which lies in the x,y=plane and is centered at the 
Origins ine surface |7]| = constant <1 is a hyperboloid of revolution 

of one sheet with an asymptotic cone whose generating line passes through 
the origin and is inclined at the angle 6 = cos”! to the Z-axis. 

The degenerate surface [yj=2 is the Z-axis. The surface 1] =0 is the 
xX,y-plane except for the circular disk %=0. The surface Y = constant 

is again the plane through the z-axis Making the angle @ with the x, 
Z-plane. The transformation from prolate to oblate coordinates is 


made by replacing a by if in the prolate equations. 


Reference for spheroidal systems: Flammer, Spheroidal Wave Functions 6 
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APPEND Dei 


SPHEROIDAL WAVE FUNCTIONS 


The particle density distribution in the vicinity of a prolate 
spheroid is describable by the prolate angle function of the first 
kind, a Cy) and the prolate radial functions of the first and 
second kinds, RO ot Coes and Rei(C,&) » respectively. 
Flammer (6| has tabulated S51’) (See Fig. II-l, which is 
taken from his data) and the radial functions Ro (¢,5) and Ro GS) 
for different values of c and the independent variables 1 and y 
The range of ° for the radial functions in his tables is, however, 
quite restricted ( § =1.0C5, 1.020, 1.0ll, and 1.077) and consequently 
it was necessary to compute further tables of the radial functions 
(Table II-1 and a 

To compute Re suc 3) and Re the following expansions 
in terms of spherical Bessel functions and Legendre Polynomials were 


obtained from Flammer's general equations. (1.1.15) for Renn (6,2) and 


(4.2.6) for Rena (6,3) 


Ror (c g) = Fa er — 5 4. PC) 4 (e$) (11-1) 
and 


7. Ol a 
Re mie oi 2 i 2. de OG, =p a At ote Cc) ae s)f (11-2) 


ye! 
where from a eae equation (1.1.8) and (l.2.5b): 


ge ($) = a v2 (9) II=la 


eee he! ~2a 
ko (c)z - = 2 dr Cc) ‘aa 


and where the prime over the summation sign indicates that the 


and 


: : O| Cc] 
summation 1s over only odd values of r, and where or (c) and dic) 
12), 





are e>pansion coefficients. Flammer tabulates these coefficients 
to five significant figures. 

The prolate radial functions of the first kind, eS (6,5) (Tabie 
II=#1), were calculated using a desk calculator and the Mathematical 
Tables Project Tables of Spherical Bessel Functions [11] . The figures 
have been spot checked, and are believed correct to plus or minus one 
im the Wace dygic. 

The Prolate Radial Functions of the second kind, K)s)(Table IT=2) 
were calculated on the CDC-160); Computer using equation ([I=2), Most 
of the Legendre Functions required for the expansions were also genera= 
ted on the computer, using the following equations given by Jahnke 


and Emde [10| . 





16 3-S-°-(2n-1) , 
r oy Net) eM? n-n-fn3) ned 
" ‘%) : a! ; “2lar-)§ * BG, Nan3) -...|(113) 


Qn (#) = Pals\Qol3) — 2 HPnel$) Rm(8) GE 


Because of the slow convergence of the equations (II-3) for large 
values of & only those functions for 2 < }.2 were computed. The 
values for é 2 1.2 were taken from Tables of Associated Legendre 
Functions [12] . The Values of Q, for $ < 1.2 were calculated on 


a desk calculator using the following relation from Jahnke and Emde [10]. 


ry 


Qolg) = cath "(Z) = 7 I1=5 
and tables of logarithms. 
The radial functions R" (63) and Roy (c,%) are tabulated here 
for values of c between 0.1 and 2.0 and for values of between 1.0 


and3.6 and 1.0 and 2.0 respectively. These functions are graphed in 


Figures (II-2), (II-3), and (II-li). The tabulated values are correct 
125 





to the number of digits shown, except that all five digit 
results are considered correct to plus or minus one unit in 

the last digit since Flammer's coefficients are listed te only 
five significant digits. The decreasing number of significant 
digits for higher values of c and > is a consequence of the 
decreasing rate of convergence of the series (II-2) for rele- 
tively large c and y e The number of significant digits was 
determined by inspection of the terms and generous estimates cf 
the maximum magnitude of the neglected terms. The author has 
endeavered to err on the conservative side in reporting these 


values, and believes all digits reported to be correct. 
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